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We discuss a topological classification of insulators and superconductors in the presence of both 
(non-spatial) discrete symmetries in the Altland-Zirnbauer classification and spatial reflection sym- 
metry in any spatial dimensions. By using the structure of bulk Dirac Hamiltonians of minimal 
matrix dimensions and explicit constructions of topological invariants, we provide the complete 
classification, which still has the same dimensional periodicities with the original Altland-Zirnbauer 
classification. The classification of reflection-symmetry-protected topological insulators and su- 
perconductors depends crucially on the way reflection symmetry operation is realized. When a 
boundary is introduced, which is reflected into itself, these non-trivial topological insulators and 
superconductors support gapless modes localized at the boundary. 



I. INTRODUCTION 

Topological insulators and superconductors arc sym- 
metry preserving fermionic systems with a bulk energy 
gap 1 ' 2 . Relevant symmetry conditions that are necessary 
to define these symmetry protected topological states can 
be divided into two categories: non-spatial symmetries 
and spatial symmetries. The Hamiltonians of non-spatial 
symmetric systems may possess time reversal symmetry 
(TRS), particle- hole symmetry (PHS), or chiral symme- 
try. They have gapless boundary states that are topo- 
logically protected and are related to physical quantities, 
such as the Hall conductivity. The subject started with 
the recognition by Kane and Mele 3 ' 4 that by incorporat- 
ing a spin-orbit coupling in the tight-binding model for 
graphene, the system will become what is now known as a 
2D Z 2 topological insulator with time reversal symmetry 
(TRS), as well as the theoretical prediction 5 and experi- 
mental observation 6 of such Z 2 topological states in the 
HgTe/CdTe quantum well. After that, three-dimensional 
(3D) Z 2 topological insulators were predicted 7-9 and 
observed; 10-13 the identification of 3 Hc-B as a topolog- 
ical superconductor was realized. 14-17 It turns out that 
those topological insulators and superconductors are just 
a part of a larger scheme — a complete classification of 
topological insulators and superconductors has been de- 
veloped with a unified periodic table 18-20 . 

Topological phases protected by these non-spatial dis- 
crete symmetries are stable against spatially homoge- 
neous as well as inhomogeneous deformations. In addi- 
tion, the protected boundary modes (edge, surface, etc) 
that appear at the boundary of topological phases in the 
periodic table are completely immune to disorder; for 
arbitrary strong disorder, as far as the bulk topological 
character is not altered in the bulk through a phase tran- 
sition, the boundary can never be Anderson- localized. 18 

With a set of discrete spatial symmetries, a topo- 
logical distinction among gapped phases {i.e., "sym- 
metry protected topological phase") can arise as well. 



One example is inversion symmetry protected topologi- 
cal insulators 21-23 where inversion symmetry is defined 
as the invariance of the system under the sign flip of the 
spatial coordinate r — »■ — r, where r = (ri,r 2 , . . . ,r<j) is 
the spatial coordinates in d spatial dimensions, r e R d 
(for lattice systems, r labels a site on a d-dimensional 
lattice). Unlike the case of non-spatial discrete symme- 
tries, for topological phases protected by a set of spatial 
symmetries, non-trivial bulk topology is not necessarily 
accompanied by a gapless boundary mode, as the bound- 
ary might break the spatial symmetries in question. The 
implications of certain specific point group symmetries 
on the topological distinction of ground states has been 
also discussed. 24-27 More importantly, Hsieh et al. 28 pre- 
dicted that one of such spatial symmetric insulators can 
be realized in the Pbi-^Sn^Te material class, and this 
prediction was experimentally verified by ARPES. 29 

In this paper, we discuss the implication of a re- 
flection (or mirror) symmetry in one spatial direction: 
it is an invariance of the system under the sign flip 
of, say, the first component of Cartesian coordinates, 
r — > r = (— n, r 2 , . . . , r^). (For an earlier study of this 
subject, see Ref. 30.) While an inversion symmetry sin- 
gles out a special point, reflection symmetry singles out 
a special (d — l)-dimensional plane (ri — in the this 
case). As a consequence, when we terminate the system 
with a (d — l)-dimcnsional boundary (plane) which is or- 
thogonal to the reflection plane {r\ — 0), the boundary 
with constant j-j (i ^ 1) is reflection symmetric under 
(ri,r 2 , . . . , rd) — > {—Vi,r2, ■ ■ ■ ,Td)- Reflection symmetry 
is arguably the simplest spatial symmetry of a system 
for which certain boundaries can respect the spatial sym- 
metry in question. This boundary property is an impor- 
tant distinction from the inversion symmetric topological 
phases, for which a plane boundary to the system alone 
does not inherit the spatial symmetry (inversion symme- 
try) in the bulk. (A pair of boundaries can be inver- 
sion symmetric to each other, though). With the special 
choice of the boundary above, we will argue that for topo- 
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logically non-trivial phases protected by reflection there 
is a stable boundary mode, in the manner similar to topo- 
logically phases protected by non-spatial discrete symme- 
tries. The correspondence still holds between the non- 
trivial bulk topology and the gapless boundary modes 
when the boundary that reflects to itself is chosen as 
shown in Figure 1. 

Although topological insulators and superconductors 
protected by non-spatial or spatial symmetries have been 
studied separately, their recognition does not directly 
provide a complete classification of topological systems 
in the presence of both non-spatial and spatial discrete 
symmetries. Therefore, we will consider the topological 
classification of reflection symmetric systems with a sub- 
set of the three non-spatial symmetries: TRS, PHS, and 
chiral symmetry. We found the topological classification 
depends not only on the set of symmetries which are 
respected but also on the way how reflection symmetry 
is realized, i.e., algebraic relations satisfied among reflec- 
tion and other non-spatial discrete symmetries when they 
exist. Our result of the classification of reflection sym- 
metric systems is summarized in Table I. The algebraic 
relations among reflection and non-spatial symmetry op- 
erations are denoted by R± and R±,± in Table I. 

Non-trivial topological states displayed in Table I are 
characterized by a topological invariant of integer (Z 2 ) 
or Z 2 type. For example, the entries in Table I marked 
by MZ indicate the presence of topologically protected 
states by a topological invariant defined on mirror in- 
variant planes in the Brillouin zone ("mirror topological 
invariant" ) . These topological states include 3D topolog- 
ical insulators protected by the "mirror Chern number" 
discussed in Ref. 30, and 2D topological superconductors 
with TRS and reflection symmetry (class Dili +R) dis- 
cussed in Ref. 31, which are characterized by reflection 
winding numbers in the ID mirror lines. We will general- 
ize those mirror numbers to any spatial dimensions and 
relevant symmetry classes. Furthermore, we will show 
that some systems are protected by the Z topological 
number and the mirror Chern number simultaneously; 
the larger one of these two numbers gives a new inte- 
gral topological invariant (denoted by Z 1 in Table I). In 
other cases, the topolgoical insulators and superconduc- 
tors in the original periodic table turn out to be invariant 
under reflection. If this is the case, the same topologi- 
cal invariant also characterizes the non-trivial topology of 
reflection-symmetric topological insulators and supercon- 
ductors. These cases are indicated by "0", Z 2 , and Z in 
Table I. In short, we claim that the reflection-symmetric 
topological states are characterized by one of the topo- 
logical invariants, '0', Z 2 , Z, MZ, and Z 1 . 

In this manuscript, we use "the minimal Dirac Hamil- 
tonian method" to characterize the AZ symmetry classes 
with reflection symmetry. Without reflection symmetry 
the topological classification (Table II) of the AZ sym- 
metry classes can be studied by Anderson localization, 18 
K-theory, 19 and minimal Dirac Hamiltonians. 32 ' 33 The 
minimal Dirac Hamiltonian method provides a direct way 



to produce the original classification (Table II). In this 
method, we first write down a bulk Dirac Hamiltonian 
preserving system's symmetries in the minimal matrix 
dimension. The topological class of the system is deter- 
mined by the existence of a symmetry preserving extra 
mass term (SPEMT), which keeps the system in the same 
topological phase during the continuous deformation. If 
this term exists in the minimal Dirac Hamiltonian, this 
phase is characterized by '0' topological invariant. If not, 
we consider a bigger system including two minimal Dirac 
Hamiltonians. The presence of a SPEMT in this system 
of the two copies implies Z 2 topological invariant charac- 
ter. Otherwise, the absence of a SPEMT implies Z char- 
acter. When classifying reflection symmetry topological 
insulators and superconductors, we study the existence 
of a SPEMT in Dirac Hamiltonians to determine topo- 
logical characters. Complementary to the minimal Dirac 
Hamiltonian method, we also look for topological invari- 
ants (0, Z 2 , Z, MZ and Z 1 ) in the presence of reflection 
symmetry to determine bulk topology. The classification 
of bulk topology in terms of topological invariants is fully 
consistent with the minimal Dirac Hamiltonian method. 

Topological insulators might have topological invari- 
ants of strong and weak indices. 34 In the original classi- 
fication table 18,19 for a d-dimcnsional system, the strong 
index is the topological invariant in d dimensions and the 
weak indices are captured by the strong indices in the di- 
mensions less than d. However, the complications arise 
when the weak indices are considered in the classification 
of reflection symmetric systems. The reason is that the 
weak indices may not be described by the strong indices 
in the dimensions less than d in the reflection classifi- 
cation table. Moreover, the weak indices might depend 
on spatial directions. That is, in different directions the 
weak indices are different because reflection symmetry 
operation only flips one direction. In this manuscript, we 
shall focus on only the strong indices of the classification. 

To name a few physically interesting topological sys- 
tems in Table I, in symmetry class All in three dimen- 
sions, with reflection symmetry specified by R , there 

are topological insulators protected by the MZ invariant. 
These are nothing but the topological insulator that was 
proposed by Hsieh et a/. 28 and observed by Xu et al. 29 
in the Pbi-^Sn^Te material class. It is the first exper- 
imental realization of crystalline topological insulators. 
This particular reflection symmetric topological insula- 
tor continues to be topologically non-trivial even in the 
absence of TRS, as indicated by "MZ" in symmetry class 
A in three dimensions in Table I. Table I also includes 
topological superconductors protected by reflection sym- 
metry, such as with TRS and reflection symmetry 2D 

topological superconductors (class Dili +R ), 31 which 

are classified by an integral- valued topological invariant. 
For symmetry class D in two dimensions, which hosts 
T-breaking topological superconductors in the absence 
of reflection symmetry, there is a reflection symmetric 
topological superconductor characterized by a Z 2 topo- 
logical invariant. Other examples will also be discussed 
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FIG. 1. For topological insulators and superconductors that 
are protected by reflection symmetry, the correspondence be- 
tween gapless surface states and bulk topology holds only 
when the surface that reflects to itself is chosen. The fig- 
ures shows that nontrivial bulk topology guarantees gapless 
states in the self-reflected surfaces. Furthermore, gapped 
states in the non-self-reflected surfaces does not imply trivial 
bulk topology. 

in Sec. VI. 

The manuscript is organized as follows: in Sec. II, 
we provide the background knowledge of reflection sym- 
metry in band theory, and in particular, describe how 
we distinguish different realizations of reflection symme- 
try operation in the presence of other non-spatial dis- 
crete symmetries. In Sec. Ill, we review the connection 
between the minimal bulk Dirac Hamiltonians and the 
topology possessing '0', Z2 and 1 topological invariants 
in the Altland-Zirnbauer(AZ) symmetry classes without 
reflection symmetry. Moreover, by considering reflection- 
symmetric Dirac Hamiltonians we show the correspon- 
dence between gapless boundary states and bulk topol- 
ogy. In Sec. IV, we consider the two kinds of the classifi- 
cations: the reflection symmetry operator commutes with 
all of the non-spatial discrete symmetries and anticom- 
mutes with TRS and PHS operators. In Sec. V, we clas- 
sify the remaining of the AZ symmetry classes possessing 
TRS and PHS under the condition that one of these sym- 
metry operators commutes with the reflection symmetry 
operator and the other anticommutes with the reflection 
symmetry operator. In Sec. VI, the concrete examples 
for topological insulators/superconductors protected by 
reflection symmetry are provided. 

II. REFLECTION SYMMETRY IN BAND 
INSULATORS 

To describe our schemes for reflection symmetric topo- 
logical insulators and superconductors, in the following 
we will start from a tight-binding Hamiltonian. We will 
focus on electronic insulators, i.e., systems with a con- 
served U(l) charge, but a similar tight-binding formalism 



can be developed for BdG Hamiltonians of topological 
superconductors. Let us consider a tight-binding Hamil- 
tonian, 

H = J2i>Hr)H(r >r ')i>(r'), (1) 

r,r' 

where ip( r ) is a iVf-component fermion annihilation op- 
erator, and index r = (n, r%, . . . , r ( i) labels a site on a d- 
dimensional lattice (the internal indices are suppressed). 
Each block in the single particle Hamiltonian H(r,r') is 
an Nf x Nf matrix, satisfying the hermiticity condition 
H^(r',r) — rl(r,r'), and we assume the total size of the 
single particle Hamiltonian is NfV x NfV, where V is 
the total number of lattice sites. The components in ip(r) 
can describe, e.g., orbitals or spin degrees of freedom, as 
well as different sites within a crystal unit cell centered 
at r. 

Provided the system has translational symmetry, 
rt(r,r') — TA.(r — r'), with periodic boundary conditions 
in each spatial direction (i.e., the system is defined on 
a torus T d ), we can perform the Fourier transformation 
and obtain in momentum space 

H=J2 ^(k)U(k)^k), (2) 

keBZ 

where the crystal momentum k runs over the first 
Brillouin zone (BZ), and the Fourier component of 
the fermion operator and the Hamiltonian are given 
by V(r) = ^ 1/2 £ feeB ze*'^(fc) and H{k) = 
J2 r e~ lk r 'H(r), respectively. The Bloch Hamiltonian 
ri{k) is diagonalized by 

n{k)\u a (k)) = e a (k)\u a (k)), a = l,...,N f , (3) 

where \u a (k)) is the a-th Bloch wavefunction with en- 
ergy e a {k). We assume that there is a finite gap at 
the Fermi level, and therefore we obtain the unique 
ground state by filling all states below the Fermi level. 
(In this manuscript, we always adjust e a (k) in such a 
way that the Fermi level is at zero energy.) We as- 
sume there are N_ (N + ) occupied (unoccupied) Bloch 
wavefunctions with N + + N- = Nf. We call the set 
of filled/unfilled Bloch wavefunctions as {|%(fc))} = 
{\va(k))}, {\u+(k))} ee {\w a (k))}, respectively, where 
hatted indices a — 1, . . . , AT_ labels the occupied bands 
only. 

In discussing symmetry protected topological phases, 
we consider a set of (discrete) symmetry conditions 
imposed on the tight-binding Hamiltonians. Altland- 
Zirnbaucr discrete symmetries, i.e., TRS, PHS, and chi- 
ral symmetry, act, on the Bloch Hamiltonian, as 

T- 1 'H(-fc)r = H(k), 
C- 1 U(-k)C = -H(k), 
S^UifyS = -H{k), (4) 

respectively, where T and C are antilinear operators, and 
S is a unitary operator. These are on-site (purely local) 
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c Spinless systems 
d Spin- 1/2 systems 

e Spin-1/2, C 2 = 1, 5(7(2) symmetry for the spin. 

f Spin-1/2, C 2 = 1 and T 2 = -1 , 5(7(2) symmetry for the spin. 



TABLE I. The complete classification table of reflection-symmetric topological insulators and superconductors: For class AIII, 
R± indicates that the reflection symmetry operator (R) commutes/anticommutes with S. For four real symmetry classes 
(R±±, -R±=f) that have TRS and PHS, the first sign ± of R indicates that R commutes/anticommutes with T and the second 
sign ± indicates that R commutes/anticommutes with C. For the four other real symmetry classes (R±) that preserve only one 
non-spatial symmetry, the sign ± indicates that R commutes/anticommutes with system's non-spatial symmetry operator. The 
Hamiltonian in the mirror symmetry plane can be block-diagonalized to two blocks in the eigenspace R = ±1. The superscript 
of 2 in the mirror symmetry classes (MSC) (See Appendix B) indicates that these two blocks of R = ±1 are independent. 



symmetries. While PHS can most naturally be intro- 
duced in the context of BdG Hamiltonians, one can still 
impose a PHS for electronic systems with conserved par- 
ticle number. 

On the other hand, reflection (1Z) is a non-local op- 
eration; by definition, a reflection 1Z in x-direction (= 
r\ -direction), say, connects fermion operators at r = 
{r-L,r 2 ,...,r d ) and at f = (~n, r 2 , . . ■ , r d ), as 

TZipi^TZ- 1 = Rifr(r), (5) 

where R is an Nf x Nf unitary matrix implementing 
reflection. The invariance of H under 1Z implies, in mo- 
mentum space, 

R- 1 H(k)R = H(k) (6) 
where k = (— hi, k 2 , ■ ■ •) = (— k\, k±). 



For example, for a spineless system, possible realiza- 
tions of these symmetries are R — I, T = 19, and 
C = t x Q 1 where t x is the first Pauli matrix acting on 
the particle-hole grading (in the BdG Hamiltonian), and 
is the complex conjugate operator. For a spin-one-half 
system, R = i<r x , 35 T = ia y Q, and C = t x Q, where Oi in- 
dicates Pauli matrices acting on spin degrees of freedom. 
We will consider more realizations of these symmetries 
later. 

Before discussing different realizations, we here note 
that, when there is a conserved U(l) charge, there is 
a phase ambiguity in the definition of the reflection 
operator; 35 when a Hamiltonian is invariant under a re- 
flection, 1Z : ip(r) -4- Rip(r), the system is also invari- 
ant under the reflection followed by a U(l) gauge trans- 
formation, ip( r ) ~^ e l ^Ri[>(f), where e 1 ^ is an arbitrary 
phase factor. The combined transformation, 1Z' : ip( r ) ~~ 
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R'ip(f) with R' = e^R, is also qualified to be called re- 
flection operation. This redefinition changes, e.g., the 
eigenvalues of the reflection transformations. 

In this manuscript, we require R to be hermitian. For 
example, in the spin-1/2 case, we add an extra 3n/2 
phase factor in R so that R — a x anticommutes with 
T and C. With this convention, we construct the classi- 
fication tables in terms of possible commutation and an- 
ticommutation relations of R with the three non-spatial 
symmetry operations. In this regard we can display the 
classification tables in a well-organized manner. While 
this is a matter of convention in classifying electrical in- 
sulators, this may not be so in classifying superconduc- 
tors (BdG systems). (It should be noted that for a AZ 
symmetry class that can be interpreted as a BdG sys- 
tems, it can also be realized as an electrical system with 
some fine tuning). 

One reason for this convention is that for different 
choices for the phase of of reflection operator (e.g., R' = 
e^R), its algebraic relation with T and C is different. 

Note that while 1Z, C, T, when acting on a fermion op- 
erator, may not commute due to the phase factor e^R, 
they always commute when acting on any fermion bi- 
linears; Physically, all of the point group symmetry op- 
erators are expected to commute with three non-spatial 
symmetry operators: TR, PH, and chiral symmetry op- 
erators. Although the requirement of the hermiticity of 
R may not correspond to the real system, it simplifies 
the classification tables. 

We will consider the topological classification when 
hermitian R commutes or anticommutes with T, C, and 
S. For simplicity, we define Rs, Rt, and Rc obeying 

SRS- 1 = R S R, TRT- 1 = R T R, CRC- 1 = R C R. (7) 

Hence, Rx = ±1 indicates the commutation or anticom- 
mutation relation between R and the non-spatial symme- 
try operator X. Furthermore, for the complex symmetry 
classes, we define the symbol of the reflection symmetry 
operator Rr s to display the algebraic relation between 
R and S. For four real symmetry classes that preserve 
TRS and PHS, the symbol Rr t .r c shows the similar 
property for T and C and provides the relation between 
R and S, which is the combination of T and C . For the 
four other real symmetry classes that have only one non- 
spatial symmetry (TRS or PHS), the symbol of the re- 
flection symmetry operator Rr t/c is defined to show the 
algebraic relation between R and the non-spatial sym- 
metry operator. In short, class AIII, AI, D, All, and 
C, which preserve only one non-spatial symmetry, have 
two possible reflection symmetry operators i?_ and R + . 
On the other hand, class BDI, Dili, CII, and CI, which 
preserve TRS and PHS, possess four possible reflection 

symmetry operators R , R++, R |_, and R + _. To 

further simplify our notations, we define in a real sym- 
metry class the reflection symmetry operator (R + and 
R++) that commutes all non-spatial symmetry operators 
as 1Z + . Similarly, 1Z~ indicates the reflection symmetry 
operator (R_ and R ) anticommuting with T and(or) 
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TABLE II. The original classification table of topological insu- 
lators and superconductors without reflection symmetry: 18,19 
The first column represents the names of the ten symmetry 
classes, associated with the presence or absence of TR, PH, 
and Chiral symmetries in the last three columns. The number 
in the last three columns denotes the absence of the sym- 
metry. The numbers +1 and —1 denote the presence of the 
symmetry and indicate the signs of the square TR operator 
and the square PH operator. 



C in a real symmetry class. 



III. TOPOLOGICAL CLASSIFICATION OF 
DIRAC HAMILTONIANS WITHOUT 
REFLECTION SYMMETRY 

To capture the essential topological features in an ef- 
ficient manner, we use the minimal Dirac Hamiltonian 
method. 33 This method simply considers the minimal 
matrix form of Dirac Hamiltonians in each AZ symmetry 
class and spatial dimension. (See below for more details.) 
When applied to topological insulators and superconduc- 
tors without spatial symmetries, this method reproduces 
the periodic table of topological insulators and supercon- 
ductors in the AZ symmetry classes (See Table II). Such a 
Dirac Hamiltonian represents a generic Hamiltonian with 
the same topological features when the spectra of the 
two systems can be continuously deformed from one to 
the other without closing the bulk band gap. Therefore, 
we will still use the method of minimal Dirac Hamiltoni- 
ans to classify topological phases of reflection symmetric 
topological insulators and superconductors. 

First, let us review the method of minimal Dirac 
Hamiltonians. The minimal Dirac Hamiltonian in d spa- 
tial dimensions which respects the set of symmetries un- 
der consideration is written as 

d 

R. = TO70 + fci7i + ^2 k Hii ( 8 ) 

where m is a constant and gamma matrices 7, obey the 
anticommutation relations 

{7i,7j} = 2<5yl. i = 0, 1, ■ ■ ■ ,d. (9) 



6 



In this manuscript, we consider the classification of 
reflection symmetric systems with three local symme- 
tries: TRS, PHS, and chiral symmetry That is, we seek 
the topological features of the Altland-Zirnbauer (AZ) 
symmetry classes 36 with reflection symmetry. When the 
Dirac Hamiltonian preserves TRS, PHS, and chiral sym- 
metry, the gamma matrices satisfy 



[fo,T\=0, {7*o,T} = 0, 
{ 7o ,C} = 0, [7*0,01=0, 

hi,s} = o, 



(10) 

(11) 

(12) 



respectively. In addition, when the system preserves re- 
flection symmetry, each term in the Dirac Hamiltonian 
obeys 



{ 7l ,i?} = 0, [7^,^=0. 



(13) 



Consider a system in a real symmetry class, which pos- 
sesses TRS or(and) PHS. If a term has the same anti- 
commutation or commutation relation of 70 in Eq. (10) 
or(and) Eq. (11) and anticommutes with the other 
gamma matrices in the Dirac Hamiltonian, we call this 
term as a mass term. Similarly, if a term behaves like 
7*0 m Eq- (10) or(and) Eq. (11), this term is identified 
as a kinetic term. Furthermore, if a mass term anticom- 
mutes with 70, this is called as an extra mass term. For 
a complex symmetry class, we also can define an extra 
mass term, which anticommutes with each gamma ma- 
trix in the Dirac Hamiltonian. 

The presence of an extra mass term in the Dirac Hamil- 
tonian plays an essential role in distinguishing topological 
phases in a system. We consider a system that preserves 
non-spatial symmetries above with/ without other sym- 
metries. If any extra mass term does not preserve the 
system's symmetries, this term cannot be added to the 
Hamiltonian. When m = 0, the bulk spectrum becomes 
gapless. With negative energies filled, we identify m = 
as a quantum phase transition point, where the gap be- 
tween the empty and occupied band closes at k = 0. That 
is, the phases with positive and negative m are topolog- 
ically different. On the contrary, an extra mass term 
preserving the system's symmetries can be added to the 
Hamiltonian as a perturbation. We define this term as 
a symmetry preserving extra mass term (SPEMT). It is 
worth to mention that by definition in a real symmetry 
class an extra mass term is always a SPEMT. When m 
varies from —00 to 00, there are no gap closing points. 
Therefore, the system for any m is always in the same 
phase, which is topologically trivial. 



A. Topological invariant '0' 

For the topological classification of the AZ symmetry 
classes there are three different kinds of topological in- 
variants: 0, Z 2 , Z. For a given set of symmetries and 
spatial dimension, we write down a Dirac Hamiltonian of 
the minimum matrix dimension, which is in the form of 



Eq. (8). If a SPEMT (9Jt) is allowed to be added to the 
Hamiltonian, the system is always in the trivial phase; we 
can classify this phase as topological invariant '0'. For 
the other two cases (Z 2 and Z) any SPEMT does not 
exist in the minimal model. 33 Therefore, the system has 
at least two different phases by varying m in Eq. (8). To 
distinguish Z 2 and Z, we need to enlarge the Hamiltonian 
and then check the presence of a SPEMT. 



B. Topological invariant 'Z2' 

While enlarging the Dirac Hamiltonian, we consider in 
the new system two minimal Dirac Hamiltonians which 
may have the same or opposite orientations. That is, one 
is given by Eq. (8) and the other is in the form of Eq. (8) 
with some ji — > —7$. Moreover, each new gamma matrix 
in the enlarged Hamiltonian must anticommute with each 
other and keep the original symmetries. The expression 
of the enlarged Hamiltonian of the two minimal Dirac 
Hamiltonians can be written as 



(14) 



remain 



The orientation of the second minimal Dirac Hamiltonian 
is determined by a z and to simplify the expression of the 
equation, let m — k . The first summation is over arbi- 
trary set of 7„ 4 (rij = 0, 1, 2, d) and the second summa- 
tion is over j n . 's that are not picked up by the first sum- 
mation. For the system with a Z 2 topological invariant, 
a SPEMT can always be added to the enlarged Hamil- 
tonian in Eq. (14) so the system is in the trivial phase. 
Therefore, the corresponding symmetries and spatial di- 
mension restrict that the system can be in the only two 
different phases when the system is characterized by the 
sign of m in the minimal Dirac Hamiltonian. 



C. Topological invariant ' 



For the system with a Z (or 2Z) topological invariant, 
when the first summation in Eq. (14) includes odd num- 
ber of 7 ni 's, a SPEMT can be treated as a perturbation 
added into the Hamiltonian. However, when there are 
even number of j ni 's in the first summation, a SPEMT 
does not exist. Therefore, the system can go through a 
quantum phase transition as m varies from positive to 
negative. 

To explain the Z(2Z) invariant, we consider n copies 
of the minimal Dirac Hamiltonian but with different m's 



7n, 



/ mi 













TO 2 



















V° 











^fc„ 3 .7i®I. (15) 



Assume all m^'s are positive. When one of m^s varies 
from positive to negative, the system goes through a 
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quantum phase transition. This phase transition can- 
not be avoided in the absence of SPEMTs due to the 
Z(2Z) invariant 33 . By adjusting m^'s the system passes 
through n times of different quantum phase transitions. 
Hence, in different raj's, n + 1 different quantum phases, 
which is labeled by the Z topological invariant, describe 
the system. 

In a system possessing a Z (or 2Z) topological invari- 
ant, we compute the topological number of the Dirac 
Hamiltonian in the form of Eq. (15). Let m, = M — k 2 . 
By calculating the winding number and the Chern num- 
ber (see appendix A), the topological number in this sys- 
tem is n(2n) if M is positive and if M is negative. As 
M is positive, one sign flipping (7, — > — 7,) causes that 
the topological number changes its sign. When the two 
systems with and without the sign switching are coupled 
a SPEMT can be present in the mixed Hamiltonian. The 
entire system in the trivial phase is consistent with the 
zero value of the topological number (n — n = 0). 



D. The correspondence between gap less edge 
states and bulk topology 

The topologically protected gapless edge states are 
present on the boundary between the trivial and non- 
trivial phases. The presence and absence of such gap- 
less edge states determines the topological bulk phases: 
(i) the presence of intact gapless edge states implies the 
non-triviality in the bulk, (ii) the absence of such states 
guarantees the system in the trivial phase. However, the 
latter statement is not always true 23 when spatial sym- 
metries are introduced. In the following, we will show 
that this statement holds when a system only preserves 
TRS, PHS, or chiral symmetry Furthermore, for re- 
flection symmetry this correspondence between bulk and 
boundaries is also true, when we choose the boundaries 
reflected to themselves and the translation symmetry in 
the reflection direction is preserved. 

First, consider the minimal Dirac Hamiltonian in 
Eq. (8) in a symmetry class with m = M — k 2 . More- 
over, we consider a domain wall in the r d direction: let 
M = M be a positive constant in the region rd > 0, 
which is in the non-trivial phase. For the trivial phase 
region r d < 0, let M — — M be a negative constant. 
Therefore. k ( j is not a good quantum number. We re- 
place kd by —id/drd- The Dirac Hamiltonian can be 
rewritten as 

d 

H=7o(mI-i7o7d^— ) + ^A;j7i, (16) 
ard i=i 

where m — M + (gf-) 2 — J2i=i ■ To have the gap- 
less energy states, we expect to find the wavefunctions 
so that the terms in the parenthesis vanish. To satisfy 
this vanishing condition, there are two possible solutions 
*7o7d0 = ±</>- We choose the minus sign to have the 



normalizable wavefunctions 

$(r d > 0) = ( Cl e-5 (1 - m - )rd + c 2 e-5 (1+m - )rd )0, 
$(r d < 0) = e-*Q- m +> T *$, (17) 

where m± = \Jl± 4Mo — 4 J2i=i kf ■ Our focus is on 

the low energy spectrum near k — so M n > X^=i &i ■ 
Therefore, m + > 1 and Rc(m_) < 1 show that the wave- 
function is normalizable. 

Because i-fold commutes with 7^0, d, by using the 
projection operator P = (I — ijojd)/2 we can discuss 
the projective Hamiltonian for the edge states in the 
(hold = -1) eigenspace 

d-l 

%sff = ^fc»7pi. ( 18 ) 
i=l 



where j pi — P7^P. The energy spectrum (±y J2i=i ^i) 
shows the gapless behavior of the edge states. Using 
the projective Hamiltonian, we can prove statement (i) 
by considering a domain wall. When both sides of the 
domain wall are in the trivial phase, a SPEMT(r) can 
be added into the Hamiltonian. Because T anticommutes 
with all of the other gamma matrix, 7 = PrP does not 
vanish and anticommutes with all of 7 P j's. Therefore, 
the gapless edge states become gapped. Furthermore, 
this statement implies that when the gapless edge states 
are intact, at least one side of the domain wall must be 
in a non-trivial phase. 

To investigate statement (ii), we focus on the behavior 
of the bulk Hamiltonian when the edge states are gapped 
without breaking any symmetry. The only one way to 
gap the edge states in V. c s is to add a symmetry pre- 
serving term that anticommutes with H c q , say 7. In the 
bulk Hamiltonian there exists corresponding symmetry 
preserving T so that 7 = PrP. Therefore, T must com- 
mutes with ijold- 37 There are two possibilities of com- 
mutation and anticommutation relations of T with each 
7j. First, T anticommutes with each 7j. 38 Second, T anti- 
commutes with 7i^o,d but commutes with 70 and jd- 39 In 
the first case T plays a role of SPEMT keeping the system 
in the trivial phase. The second case needs to be investi- 
gated scrupulously. Although T is not a SPEMT, in the 
case of some specific symmetries there exist a SPEMT, 
which is a hermitian matrix iTjo-jd- If the system pre- 
serves TRS, PHS, and chiral symmetry, then by Eqs. (10) 
to (12) iT-fQ-fd also preserves those symmetries. The cor- 
respondence between edge states and bulk topology can 
be applied for these three symmetries. 

Let reflection symmetry reflect only in the kd direction, 
then {i?,7d} = and [R,^ d ] = 0. Therefore, if^ d 
breaks reflection symmetry and then the absence of the 
gapless edge states does not imply the trivial bulk topol- 
ogy. However, if reflection symmetry is not in the k d 
direction, then iT"f jd preserves the symmetry and can 
be present in the Hamiltonian as a SPEMT. The gapped 
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edge states possessing reflection symmetry guarantee the 
triviality in bulk. 

When the translational symmetry in the reflection di- 
rection is broken, the correspondence between gapless 
edge states and bulk topology does not hold. However, 
we still can use the mid-gap states in the entanglement 
spectrum to distinguish topological trivial and non-trivial 
phase. 21,23 We leave this issue in the future discussion. 
In the manuscript, we use the existence of SPEMTs in 
the minimal Dirac Hamiltonians to determine possible 
topological phases. 

IV. THE CLASSIFICATION OF R+, K + AND TZ~ 
- SHIFTED PERIODIC TABLE 

We consider the real AZ symmetry classes with the re- 
flection symmetry commuting (1Z + ) and anticommuting 
with T and C. For 1Z + , the classification table is 
obtained from the original table without reflection sym- 
metry by "upward shift" in spatial dimensions as shown 
in Table III. The topological invariant Z is replaced by 
a new topological invariant MZ, which will be explained 
later ("mirror" topological invariant). Similarly, for 1Z~ , 
in d dimensions the topological invariants in the new ta- 
ble is the ones in d + 1 dimensions in the original table 
( "downward shift" ) , except for the absence of the second 
descendant Z2 40 of Z as shown in Table IV. In the fol- 
lowing, for the real symmetry classes, we will construct 
these two tables for 1Z~ and 1Z + and define the topologi- 
cal invariant MZ. We leave the discussion of the complex 
symmetry classes (R+) for interested readers. 



A. Classification of Dirac Hamiltonians 

Let us start by giving a brief description of the mech- 
anism behind these dimensional shifts. By knowing that 
7! anticommutes with all of the other gamma matrices, 
we construct a hermitian matrix 171 R satisfying the an- 
ticommutation relation 

{ij 1 R,'H} = 0. (19) 

First, for the case of 11', from Eqs. (10), (11) and (13) 
ijiR can be used as a gamma matrix to construct an- 
other Dirac kinetic term in one higher dimensions (d+l), 
say as 7<j+i- Because [ijiR, C] = and {i-fi R,T} = 0, 
the (d + l)-dimensional Hamiltonian preserves the same 
set of local AZ symmetries. Alternatively, to construct 
a Hamiltonian in d dimensions with 1Z~ , we can start 
from a system in d + 1 spatial dimensions preserving the 
same local symmetries, but not reflection. By removing 
one gamma matrix 7^+1 (and momentum component) 
from the kinetic, we obtain the <i-dimensional Hamilto- 
nian with reflection symmetry. We will later make use 
of the topological classification of the d + 1-dimensional 
Hamiltonians without reflection to discuss the topologi- 



cal classification of our d-dimensional target Hamiltoni- 
ans with reflection symmetry. 

On the other hand, for the case of 1Z + , ijiR can be 
used as an extra mass term: it can be added to the 
Hamiltonian without changing its AZ symmetry class 
since {ijiR, C} — and [ijiR, T] — 0, while 171 R breaks 
the reflection symmetry. Because of the algebraic struc- 
ture of the Clifford algebra, adding a mass term effec- 
tively acts as removing one kinetic gamma matrix, and 
therefore, effectively decreases the spatial dimension by 
one. 33 Therefore, the topological classification of the d- 
dimensional Hamiltonians with 1Z + is related to the clas- 
sification of (d— l)-dimcnsional Hamiltonians in the cor- 
responding AZ symmetry class without reflection symme- 
try. This "upward" shift is also supported by considering 
the Hamiltonian in the (d — l)-dimensional mirror plane 
in the Brillouin zone. The topological invariant defined 
for the (d — l)-dimensional Hamiltonian directly deter- 
mines the topological class of the original d-dimensional 
Hamiltonian with reflection symmetry. (See below for 
more details). 

For those two cases (n~ and TZ+), when a SPEMT 
in d ± 1 dimensions in the AZ symmetry class without 
reflection symmetry is still a SPEMT in d dimensions 
with reflection symmetry, both of the systems share the 
same topological invariant. In the following we will show 
that when systems in d ± 1 dimensions possesses '0' and 
Z2 topological invariants, the corresponding reflectional 
systems have the same topological invariants. Likewise, 
for a Z invariant in d ± 1 dimensions, the corresponding 
reflection symmetric system in d dimension has a Z-like 
topological invariant. 

a. Topological invariant '0' Consider a system in an 
AZ symmetry class in d± 1 that has a '0' topological in- 
variant. Therefore, the presence of a SPEMT (9Jt) in the 
minimal Dirac Hamiltonian in Eq. (8) keeps the system 
in the trivial phase. We define the reflection operator 
R = IA71, where A = 71(7^+1) corresponds to 1Z + (1Z~). 
Because 371 anticommutes with A and 71 , the system in d 
spatial dimensions with 9Jt preserves reflection symmetry. 
The same AZ symmetry class with reflection symmetry 
in d dimensions has a '0' topological invariant. 

b. Topological invariant "L2 ' For a Z2 topological 
invariant, in d± 1 dimensions the minimal Dirac Hamil- 
tonian in Eq. (8) in one of the AZ symmetry classes has 
no SPEMTs. Therefore, SPEMTs do not exist for the 
minimal Dirac Hamiltonian in d dimensions with reflec- 
tion symmetry R = iA-fi. Again to find the topological 
property for the reflection symmetry, the minimal Dirac 
Hamiltonian in d dimensions can be enlarged in several 
ways: 

%i = fci7i ® I + ^2 ® <*z + J2 kn i^ n i I ' 

rii^l remain 

(20) 

rii^l remain 

(21) 
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with the unchanged reflection symmetry operator R = 
iAji ®I. Because of the Z 2 topological invariant, for "H 2 
we can construct the Hamiltonian in d ± 1 dimensions in 
the same AZ symmetry class with a SPEMT 971 without 
the reflection symmetry as 

n^ 1 = ni + \k®i + m, (22) 

where A = kd+i(m) in d+ l(d— 1) dimensions. To check 
whether 97t preserves the reflection symmetry, the com- 
mutation relation between R = iA-fi ® I and 97T should 
be considered. By the definition of a SPEMT, 971 anti- 
commutcs with A <£> I and 71 <8 I. Therefore, 971 preserves 
the reflection symmetry. The Hamiltonian can be gapped 
out without breaking any symmetry. 

For the second Hamiltonian "H 2 i the corresponding 
gapped Hamiltonian in d ± 1 dimensions is written as 

Ui ±l ' = Ut + AA ® a z + 971. (23) 

Because the SPEMT 97t commutes with 71 <g> I and A ® I, 
971 preserves the reflection symmetry. This is so since if 
one of 71 (8)1 and A® I anticommutes with 97t, I<8cr z must 
commute with 971 so both of 7i<8l and A (8)1 anticommute 
with 97T. This contradicts with the assumption that there 
are no SPEMTs in the minimal Hamiltonian. In short, 
the d-dimcnsional system with reflection symmetry has 
a Z 2 invariant inherited from the d ± 1-dimcnsional sys- 
tem without reflection symmetry. However, there is an 
exception. For the anticommutation case a reflection sys- 
tem from d + 1 dimensions corresponding to the second 
descendants 40 Z 2 of Z has '0' topological invariant in- 
stead of Z, which will be discussed later. 

c. Topological invariant *Z ' Consider a system in 
d ± 1 dimensions in an AZ symmetry class that has a 
Z topological invariant. Therefore, in the correspond- 
ing d-dimensional system with reflection symmetry any 
SPEMT does not exist for the minimal Dirac Hamilto- 
nian, which is similar with the Z 2 case. To distinguish 
the topological invariant from Z 2 , we need to enlarge 
the minimal Dirac Hamiltonian in the forms of Eqs. (20) 
and (21). The corresponding Hamiltonians without re- 
flection symmetry are written in the forms of Eqs. (22) 
and (23) respectively so that a SPEMT 971 preserves the 
reflection symmetry. However, by the definition of the 
Z topological invariant, the SPEMT 971 are present in 
Eqs. (22) and (23) only when the first summation is over 
odd number of the gamma matrices. Thus, for the sys- 
tem with reflection symmetry the presence/absence of a 
SPEMT is determined by the first summation odd/even 
number of gamma matrices in Eqs. (20) and (21) but 
does not depend on the way to enlarge 71. In the next 
paragraph we will prove that such a system possesses a 
topological invariant. We label this invariant by 'MZ'. 
The reason is that the MZ system behaves the same with 
the Z one in the non-reflectional symmetry direction but 
topological property of MZ is insensitive in the reflec- 
tional symmetry direction. Moreover, we will show that 
MZ number is defined in the refiection(mirror) symmetry 
planes (fci = or it) so M means mirror. 



Commutation relations of R 
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TABLE III. The classification table for i?(class A), R + (class 
AIII), and lZ + (re&\ symmetry classes). Each non-spatial sym- 
metry operator commutes with R. For class A, a system has 
only reflection symmetry so no commutation issue is in this 
class. If we treat MZ(2MZ) as Z(2Z), this table is obtained 
from the original table just by "upward shift" in spatial di- 
mensions. 



Anticommutation relations of R 



AZ class+R\rf 


1 


2 


3 


4 


5 


6 


T 


c s 


AI 








2MZ 








Z 2 


+ 





BDI 











2MZ 








+ 


+ 1 


D 


MZ 











2MZ 








+ 


Dili 


z 2 


MZ 











2MZ 




+ 1 


All 





Z 2 


MZ 
















CII 








z 2 


MZ 










- 1 


c 


2MZ 








z 2 


MZ 








- 


CI 





2MZ 








Z 2 


MZ 


+ 


- 1 



TABLE IV. The classification table for 1l~ . For the eight 
real symmetry classes, R anticommutes with T and C but 
commutes with S. For class AIII, the only non-spatial sym- 
metry operator S anticommutes with R. We treat MZ(2MZ) 
as Z(2Z), then the topological invariants in this table in d di- 
mensions is the ones in the original table in d+ 1 dimensions, 
except for the second descendant Z 2 of Z. 



B. Topological numbers 

1. Topological numbers of MZ 

For the case where the reflection operator R commutes 
with the non-spatial discrete symmetries, the topological 
numbers MZ can be defined from the bulk Hamiltonian 
to characterize bulk topology and protected gapless edge 
states. The Hamiltonian without k\ commutes with R; 
therefore, the Hamiltonian can be block-diagonalized in 
the two eigenspaces R — ±1 because R is hermitian. 
Each individual block Hamiltonian is (not) invariant un- 
der the original symmetries if such symmetry operators 
(anti)commute with R. However, the Hamiltonian still 
belongs to one of the AZ symmetry classes, which corre- 
sponds to the non-spatial symmetry operators commut- 
ing with R. We name this symmetry class in the mirror 
planes of k\ = or 7r as a mirror symmetry class. The 
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details of mirror symmetry classes will be discussed in 
appendix B. 

We focus on one of the blocks with a definite eigen- 
value, R = 1, say, because the topological numbers of 
these two blocks differ by signs when the weak topolog- 
ical index vanishes. The reason is that a d-dimensional 
system with a nonzero weak index can be understood by 
a stacking limit of d — 1 dimensional topologically non- 
trivial layers . Furthermore, for any k\ the sum of these 
two topological numbers is invariant; hence, if this total 
number does not vanish, by definition the weak index is 
nonzero. In the following, we always consider the case 
that the weak index vanishes to define the MZ number. 

a. MZ for 1Z + First, suppose the non-spatial sym- 
metry operators commute with R. For a d-dimensional 
system possessing a MZ topological invariant, the topo- 
logical number does not depend on the k\ direction (di- 
rection of the reflection symmetry). Furthermore, in d — 1 
dimensions the mirror symmetry class, which is the same 
with the original AZ symmetry class, has a Z topological 
invariant. Hence, to obtain the MZ number, we can cal- 
culate the Z number in one of the blocks with a definite 
reflection eigenvalue in d— 1 dimensions by using Eq. (A3) 
or Eq. (A6) because the Dirac Hamiltonian without 71 
commutes with R. This Z property is protected by the 
corresponding block diagonal non-spatial symmetry op- 
erators since R commutes with these operators. 

In the continuum model the MZ number can be prop- 
erly defined for the block Hamiltonians at k\ — 0. How- 
ever, in the lattice model, which can be obtained by the 
replacement k\ — > smnk\ (n € Z), sinnfci vanishes in the 
Hamiltonian only when k\ = 0, ±w/n, ±2ir/n, . . . , ±tt. 
These points are the possible positions to have a (d — 1)- 
dimensional Z number. However, the Z numbers which 
are not at the symmetry points (fci =0, 7r), are fragile, or 
not protected: They can vanish by coupling the opposite 
Z numbers in the other block of R = —1 without break- 
ing reflection symmetry. Furthermore, no SPEMTs are 
allowed in the bulk Dirac Hamiltonian around the sym- 
metry points so the Z numbers at the symmetry points 
are invariant. Therefore, we can calculate the two num- 
bers z^q" 1 and v^~ l at ki = 0, 7r respectively in the block 
ofi?=l. 

To have the topological number of the strong index, we 
need to consider translational symmetry breaking. The 
presence of translational symmetry breaking along the r\ 
direction connects the Z numbers at the two symmetric 
points so the total topological invariant number is vf^ 1 + 
v^~ x . However, this number is not the strong index. The 
strong index Nmz is given by 

N M z = vt 1 + 4' 1 - 2Ar wcak, (24) 

where iV woak is the mirror weak index, which is the weak 
index in one of the blocks of R = ±1 and invariant for 
any k\ . Such a mirror weak index is determined by d — 1- 
dimensional nontrivial layers 41 , which are stacked to a 
nontrivial weak system. To have the strong index, we 
determine the mirror weak index first by considering two 



possible situations: x v^ 1 > and 1 i>% 1 < 0. On 
the one hand (vq~ 1 v^ 1 > 0), the mirror weak index is 

iVwcak = sgn(i/ d - 1 )min(|i/ d - 1 |, I^Q. (25) 

Because the total invariant number are the sum of 2iV wea k 
and the strong index, we can write the strong index topo- 
logical number is 

sgnK^K" 1 -^- 1 !. (26) 

On the other hand, when v^~ 1 v^~ 1 < the mirror weak 
index is absent. The strong index is the total invariant 
number v^ 1 + v^~ x . From these two cases, the MZ 
number is defined as 

N MZ = sgn(^- 1 - vt'KWt'l li^l). (27) 

The signs determine the orientation; however, the Nmz 
does not have the summation property like Z topologi- 
cal invariant (Nz). Consider a system is the collection of 
several subsystems. Each subsystem has its own Z num- 
ber N z . Therefore, the Z number of the entire system 
is given by Nz = J2 t N%. This relations does not hold 
for MZ. To obtain Nmz, we have to compute the two 
Z numbers (vq -1 and z^ _1 ) for the entire system in the 
reflection symmetric planes and then use Eq. (27). 

b. MZ for TlT Similarly, we consider the case that 
a d-dimensional system with reflection symmetry opera- 
tor 1Z~ , which anticommutes with the TRS operator or 
PHS operator. If d is even, the system possessing MZ 
preserves Chiral symmetry. This Chiral symmetry op- 
erator commutes with TlT so in the Hamiltonian block- 
diagonalized by 1Z~ each block has the corresponding 
Chiral symmetry operator. Neither the TRS nor PHS is 
block-diagonalized at the same time because of the an- 
ticommutation relations. Therefore, each block in the 
Hamiltonian belongs to class AIII in even dimensions 
and then the topological number can be evaluated as the 
winding number by Eq. (A3). On the other band, if d 
is odd, the off-diagonal Hamiltonians, which do not pre- 
serve any symmetry, belong to class A in odd dimensions. 
Hence, in class A the Chern number in Eq. (A6) can char- 
acterize the topological number at the symmetry points. 
In short, the topological number in the anticommutation 
case still can be described by Eq. (27). 



2. Topological numbers 0/Z2 

a. Z2 for 1Z + Consider a <i-dimcnsional Z2 system 
with 1Z + commuting with the local symmetries. In such 
a system, the topological invariant for a block in the R- 
block-diagonal Hamiltonian is Z2 in d — 1 dimensions 
in the mirror symmetry class, which is the same with 
system's symmetry class. Such a Z 2 topological invariant 
was already evaluated in several ways 3 ' 17,34 ' 42 . Therefore, 
at the two symmetry points ki — 0, ir, the Z2 numbers 
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are defined as vq and respectively in the block of the 
Hamiltonian. The Z 2 number for the entire system is 

Nz 2 =v + Vk mod 2. (28) 

The reason is that the reflection symmetry does not 
prevent a translational symmetry breaking density wave 
from coupling and gapping out a pair of bulk band gap 
closing (quantum phase transitions) at these two sym- 
metry points. Only one bulk band gap closing survives 
under arbitrary symmetry preserving perturbations when 
a system possesses an odd number of closing. 

b. Z 2 for 1Z~ Consider in the case that reflection 
symmetry operator anticommutes with T and(or) C, 
which is more complicated. We discuss two possible cases 
respectively: the first and second descendants Z 2 of Z in 
d + 1 dimensions. 

First, consider a system with reflection symmetry cor- 
responding to the first descendant Z 2 of Z in d + 1 di- 
mensions. We note that in d dimensions the original 
topological classification gives a Z 2 topological invariant. 
The topological number in this case can be defined by the 
original Z2 number. That is, a system in such a symme- 
try class with and without the reflection symmetry has 
the same topological invariant. 

Secondly, in a system with reflection symmetry corre- 
sponding to the second descendant Z 2 of Z in d+1 dimen- 
sions, the topological number cannot be properly defined. 
The reason is that T and(or) C anticommute with R so 
the mirror symmetry class is class A. Therefore, no Z 2 
topological numbers can be defined at k\ = 0, ir. Fur- 
thermore, in d dimensions the corresponding topological 
invariant is '0' in the original classification. It turns out 
that a SPEMT can be present in the Hamiltonian to pre- 
vent the bulk gap closing so this case is classified as '0'. 
The further discussion is in the following. Without en- 
larging the minimal Dirac Hamiltonian in d dimensions 
in the corresponding symmetry class, d+3 kinetic gamma 
matrices (71,72, • • ■ ,7^+3) and one mass gamma matrix 
(70) can be present. 33,43 The reflection symmetry opera- 
tor is defined as R = i7o7<2+i to satisfy the anticommu- 
tation relations with TRS and PHS operators. Due to 
the presence of 7^+1, jd+2, ld+3, we have more choices 
to add some symmetry preserving terms in the minimal 
Dirac Hamiltonian in Eq. (8) 

d 

H s = m 7o + fci 7 i + hli + SA, (29) 

where 5 is a positive constant and A = i7i7d+i7d+2, 
which is invariant under all system's symmetries. Since 
A commutes with only 71 in H5, the eigenvalues of 
71 (fcil + 5i"id+i"id+2) arc k\ ± 5 and — k\ ± 5 due to the 
eigenvalues ±1 of ijd+iJd+2- Therefore, when the quan- 
tum phase transition (m = 0) occurs, the bulk gap closing 
points shift k\ = ±5 and k± = 0. Now we can add an- 
other symmetry preserving term to prevent the bulk gap 
closing at the new transition points by breaking trans- 
lational symmetry. This gap opening term is written in 



the form of the second quantization 

* = E Hr+n+sKck^+s + h.c.) (30) 

— ?7<fel <7J 

+ E 

— r;<fei<r; 

where 77 is a positive constant less than S and 9t = 
*7<i+i7(i+27d+3, which anticommutes all of the terms in 
Us. Also, 9T preserves TSR and PHS by Eqs. (10) 
and (11). By the definition of the reflection symmetry 
operator R — ^2 kl c\. («7o7d+i)c-fe 1 44 , it is easy to check 
that 91 preserves the reflection symmetry. The last thing 
we need to verify is that prevents the bulk gap clos- 
ing. To have the low energy spectrum, consider that 
case m = and kj_ — so Us is a function of k\. The 
Hamiltonian %s{k\) with c9T is in the form of the second 
quantization is written as 

W«,C= E (*fc 1+ «^(fcl 

— 7j<fei<r; 

+ high energy terms, (31) 
where = (c Pl+n c Pl - v ) T and 

We compute the eigenvalues of the two blocks ^H v {k\+S) 
and W^—ki —5)) to capture the low energy spectrum, be- 
cause those two blocks are the reflection symmetry part- 
ners sharing the same energy spectrum. Therefore, con- 
sider the energy spectrum of one of the blocks, say 

Hr,(ki +5) = (fci + 5)1 ® 71 + r]cj z ® 71 

+ SI <g> A + co- y ® 01. (33) 

We note that *Tt anticommutes with A and 71 and A com- 
mutes with 71. Therefore, the expression of the energy 
square is 

E 2 = (ri±yjkl+c?) 2 , (77 ± ^(^1+ 2(5) 2 +c 2 ) 2 . (34) 

Hence, when c is larger than 77, the energy never be- 
comes zero. The bulk gap closing has been blocked by 
the symmetry preserving terms as SPEMTs. Therefore, 
the topological invariant in this case is '0'. In sec. VIC, 
one example is provided to show that the translational 
symmetry breaking 91 gaps the edge states and destroys 
the mid-gap states in the entanglement spectrum. 

V. THE CLASSIFICATION OF R_+ AND 

In this section, we consider the topological classifica- 
tion of insulators and superconductors for the cases of 
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R |_ and To have these anticommutation and com- 
mutation relations we have to consider the AZ symmetry 
classes that preserve both TRS and PHS — class Dili, 
CII, CI, and BDI. The results are summarized in Tables V 
and VI. 



A. Classification of Dirac Hamiltonians with 

R = i^S 

As a start, let us consider, as a possible defini- 
tion of reflection operator, R = i"/iS. The commuta- 
tion/ anticommutation relations of R — i^\S with TRS 
and PHS operators are summarized in in Table V. They 
can be verified as follows. Let us go back to the expres- 
sions of TRS and PHS operators 

T = U T Q, C = U C Q, (35) 

where Ut and Uc are complex matrices. To simplify our 
problem, we assume that Ut and Uc are hermitian and 
unitary. To define chiral operator S, which is hermitian, 
we let S = TC if [U c , U T ] = or S = iTC if {U c , U T } = 
0. Therefore, R = ij\S is hermitian. To determine the 
commutation and anticommutation relations of R with 
T and C, we have to check the relations of S with T and 
C. In the both cases ([U c , U T ] = and {U c , U T } = 0), 
we have the same relations 

TST- 1 = ±S, CSC' 1 = ±S, (36) 

where we pick up the plus sign in front of S when T 2 = 
±1 and C 2 = ±1, whereas we pick up the minus sign 
when T 2 = ±1 and C 2 = =Fl. The reason is that 

±1 = T 2 = U T U T , ±l = C 2 = U C U* C . (37) 

By using hermitian and unitary properties of Ut and Uc, 

U T = ±U T , U c = ±U* C . (38) 

By Eqs. (10) and (11), we obtain the relations exactly 
shown in Table V: [T, R] = and {C, R} = when T 2 = 
±1 and C 2 = ±1 and {T,R} = and [C,R] = when 
T 2 = ±1 and C 2 = Tl- 

This construction of R docs not require any new 
gamma matrix. Hence, the absence of a SPEMT in 
the original classification is unchanged when the reflec- 
tion symmetry is considered. When a SPEMT exists in 
the original Hamiltonian, we have to check whether this 
SPEMT preserves the reflection symmetry. If so, the 
system is in the trivial phase. If not, the reflection sym- 
metry provides new topological phases. We consider the 
following three cases separately. 

• Let us first suppose a system in an AZ symmetry 
class without reflection symmetry has '0' topologi- 
cal invariant. A SPEMT (9Jt) exists in the minimal 
Hamiltonian. That is, 371 anticommutcs with S and 



Class 


TR 


PH 


Ch 


R 


d = 2 


d = 3 


d = 4 


d = 5 


d = 6 


AIII 








1 


R- 





Z 1 





Z 1 





BDI 


+1 


+1 


1 


R+- 











2Z 1 





Dili 


-1 


+1 


1 


R-+ 


Z 2 


z 1 











CII 


-1 


-1 


1 


R+- 





z 2 


Z 2 


z 1 





CI 


+1 


-1 


1 


R-+ 





2Z 1 





z 2 


z 2 



TABLE V. The classification table for the case of reflection 
symmetry operator given by i^/iS. By using the similar dis- 
cussion of R — iy-iS, a reflection symmetric system in class 
AIII possesses a Z invariant in odd dimensions and '0' 

invariant in even dimensions. The AZ symmetry classes with 
the reflection symmetry have Z 2 and Z 1 corresponding to Z 2 
and Z without the reflection symmetry respectively. 



7j. Therefore, when the reflection symmetry is im- 
posed on the system, R = i'jiS commutes with 371 
so the reflection symmetry is preserved. The sys- 
tem is still in the trivial phase and classified by the 
'0' topological invariant. 

• For the AZ symmetry class that has a Z 2 topo- 
logical invariant, no SPEMTs are allowed in the 
minimal Hamiltonian. However, when the size of 
the Hamiltonian is doubled in the form of Eq. (14), 
without considering reflection symmetry a SPEMT 
(SDt) does exist. Wc found that j ni ® I in Eq. (14) 
must commutes with SOT because if not, the min- 
imal Hamiltonian can have a SPEMT. Also, we 
know {-f n . eg) I, £01} = 0. Therefore, if m ^ 1, 
[R = i'jiS, OJt] = shows that the reflection sym- 
metry is preserved. Otherwise, the reflection sym- 
metry is broken. Hence, any SPEMT seems to be 
absent in this case; however, later we will show that 
a SPEMT can be present, which is similar to the 
counterfeit 1^ case in sec. IV B 2. Hence, the reflec- 
tion symmetry does not provide any extra topolog- 
ical phase. Such a system with the reflection sym- 
metry is still classified as 

• Finally, a system with a Z topological invariant 
guarantees no SPEMTs in the minimal Hamilto- 
nian. After doubling the size of the minimal Hamil- 
tonian, the original Z topological invariant forbids 
any SPEMT in Eq. (14) when there are even terms 
in the first summation. When reflection symme- 
try is considered, by the similar reasoning as the 
previous discussion, {R = i-fi S, 371} = breaks re- 
flection symmetry as = 1 in Eq. (14). In short, 
the system can be gapped out by a SPEMT only 
when the first summation in Eq. (14) is over an 
odd number of the gamma matrices, excluding 71. 
Let us go back to the discussion of the MZ topo- 
logical invariant. Eqs. (20) and (21) provide all of 
the possible twice-as-big minimal Dirac Hamilto- 
nian. Only when the first summation of Eq. (20) 
is over an odd number of the gamma matrices, a 
SPEMT can be present in the Hamiltonian. There- 
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fore, for an even number of the gamma matrices in 
the first summation of Eq. (20) a SPMET is forbid- 
den by the Z topological invariant. In Eq. (21) since 
fci7i <S> <J Z in the Hamiltonian with an even number 
of gamma matrices in the first summation, a mass 
term preserving the non-spatial symmetries breaks 
the reflection symmetry. Therefore, the bulk topol- 
ogy is non-trivial for both of the (the even numbers 
of gamma matrices). For the former, the Z number 
is non-zero as well as for the latter, the MZ number 
does not vanishes. Thus, the bulk topology might 
be protected by Z and MZ topological invariants. 
In the next subsection, we will show that for such 
a system the topological invariant labeled by Z 1 is 
determined by Z and MZ topological numbers. 

The above considerations lead to the classification 
summarized in Table V. In the following, we present 
a proper definition of topological invariants. 



1. Topological number of Z 1 

To define the proper Z 1 number, characterizing bulk 
topology and intact gapless edge modes, we consider Z 
number (JVz) and MZ number (Nmz) as our candidates. 
In the following, we will prove that the Z 1 number is 



N z i =Max(|iVz|,|JVMz|). 



(39) 



To simplify the problem, we consider the presence of 
non-trivial topology only at k\ = plane. In other words, 
Z and MZ numbers can be determined by the Hamilto- 
nian around this symmetry plane. We leave the gen- 
eral proof in appendix C for interested readers. In a d- 
dimensional system the Z number Nz and the MZ num- 
ber Nmzi which is the Z number Vd-i at &i = 0> are both 
defined by either Eq. (A3) or Eq. (A6). We define (±,±) 
and (±, =p) to describe the signs of the gamma matrices 
in Eq. (8). The first slot of the parentheses indicates 
the sign of 71 and the second slot +(— ) shows even(odd) 
number of the other gamma matrices having the minus 
sign. According to the Z 1 properties, the system with 
(+, +) and (+, -) can be gapped out by a SPEMT. Sim- 
ilarly, the pair of ( — ,—), (— , +) are in the trivial phase, 
too. However, two minimal Hamiltonians from these two 
different pairs are protected by the Z 1 topological invari- 
ant. We can treat these two pairs be two independent 
systems of (+, ±) and (— ,±). Let the numbers of the 
minimal Hamiltonians of (±, ±) and (±, =f) be iV± ; ± and 
N±^ respectively. A nontrivial minimal Hamiltonian 
provides one protected gapless edge mode so the num- 
bers of intact gapless edge modes for these systems are 
given by respectively 



N+,± = N+,+ - iV+,_, A_ 



jV_ + -AT__. (40) 



m = M — k 2 in Eq. (8). Therefore, the numbers of the 
intact gapless modes in the expression of N% and Nmz 
are written as 



N. 



N z + N, 



+,± 



Nz — Nmz 
2 



(41) 



We note that a system (+, ±) has N% — ±1 and Nmz = 
±1 and with (— , ±) has N% = =pl and Nmz — ±1 when 



The total number of the intact gapless edge modes in the 
system is the Z 1 number A z i = |AT +i ±| + |AT_.±| as in 
shown Eq. (39). 



2. Topological number of Z2 

In Table V the Z2 number from a system with reflection 
symmetry can be computed in the same way of the Z 2 
number without reflection symmetry. The reason is that 
the reflection symmetry does not give rise to any new 
topological phase. From the previous discussion, only for 
the two-copy minimal Dirac Hamiltonian "H 2 in the form 
of Eq. (21) the bulk gap closing cannot be prevented by 
any symmetry preserving homogeneous term. However, 
we will show later that this Hamiltonian can be kept 
gapped by some inhomogeneous SPEMTs. 

This Z 2 class can be separated into two slightly dif- 
ferent cases represented by the original Z 2 invariants in 
table II: the second descendant Z 2 of Z in odd dimensions 
and the first descendant Z 2 of Z in even dimensions. In 
the following, we consider both of the cases together, in 
the twice-as-big minimal Dirac Hamiltonian there exists 
a SPEMT preserving reflection symmetry (R = ijiS). 
Therefore, this reflection symmetry still keeps the origi- 
nal Z 2 invairants. 

For the second descendant, without enlarging the di- 
mension of the minimal Dirac Hamiltonian d + 2 kinetic 
gamma matrices 71 , 72, ... , 7,^+2 and one mass matrix 
7o 33 can be used for the SPEMT construction. 

For the first descendent in even dimensions, we have 
the same gamma matrices, except for missing jd+2- 
Therefore, from those gammas matrices some symmetry 
preserving terms can be constructed and added into the 
enlarged minimal Dirac Hamiltonian 

Ui = fci7i <g> a z + 6A + ^ k n % ln t ® CT 2 

+ *nj7n, ®I, (42) 



where A = (2)7^+1 Iln^™i 7"; ® a u when the first sum- 
mation is over an even number of terms or A = 
W n«j^i ln,i®o~u when it is over an odd number of terms. 
The presence or absence of i keeps A being hermitian and 
choosing u = x,y lets A preserve TRS and PHS. We can 
check that A preserves all of the system's symmetries, 
including reflection symmetry R = i-yi S ® I. The situ- 
ation is similar with Eq. (29): the presence of A shifts 
the bulk gap closing points at m = 0, k\ = ±8, k^_ = 0. 
However, this system of the twice-as-big minimal Dirac 



14 





1 it 


P 77 
111 




T> 
It 


A o 

Q — Z 


A Q 

a — o 


(A — 4 


A C 

a — o 


A f: 

a — o 


A 7 


BDI 


+1 


+1 


1 


R-+ 





2MZ 





2Z 





2MZ 


Dili 


-1 


+1 


1 


R+- 





2Z 





2MZ 





2Z 


CII 


-1 


-1 


1 


R-+ 





2MZ 





2Z 





2MZ 


CI 


+1 


-1 


1 


R+- 





2Z 





2MZ 





2Z 



TABLE VI. The classification table for the rest of the commu- 
tation and anticommutation relations. Non-trivial topology 
only shows up in odd spatial dimensions. 



Hamiltonian still can be gapped out by a special SPEMT 

— ?7<fei <7J 

+ X! (icki+v-sMc^-^-s + h.c), (43) 

— 77 < fe i < 

where 9? = (2)71 IIn^™i Tni^o'u when the first summation 
in Eq. (42) is over an even number of terms, or 91 = 

(*)7d+i7i Ilrii/i 7«i ® a u when the summation is over 
an odd number of terms. Therefore, 9t anticommutes 
with all of the terms in Eq. (42). Again the presence 
or absence of i guarantees hermiticity of 9T. By properly 
choosing u = x, y, 9? preserves TSR and PHS by Eqs. (10) 
and (11). By the definition of the reflection symmetry 
operator R = ^2 kl c\. i (i-fiS ® I)c_fe 17 91 preserves the 
reflection symmetry. Following the similar discussion in 
sec. IV B 2, the low energy spectrum is shown in Eq. (34). 
When c > i], the original bulk gap closing is blocked. 
No quantum phase transitions, no topological non-trivial 
phases. The system of the twice-as-big minimal Dirac 
Hamiltonian is trivial. Hence, this case is classified as 
Z 2 from the original topological invariant without the 
reflection symmetry. 



With extra mass term, R can be constructed and sat- 
isfy the required commutation and anticommutation re- 
lations: 

d 

R + _=ij j j 7i, asd = 4n + 2, (44) 

i=0,»^l 
d 

R-+ = 7 n 7», as d = 4 "- (45) 

i=0,i#l 

It is easy to check that 7 preserves the reflection symme- 
try. Hence, this gamma matrix plays a role in a SPEMT. 
This shows that as R T = R c = ±1, d = 4n(+2), 
those four symmetry classes possess the '0' topological 
invariant, which explains several 0's in even dimensions 
in Tables V and VI. 

To have the opposite commutation and anticommuta- 
tion relations in the previous case, we use "fd+i to con- 
struct R 



d+l 

r-+=i n 

i=0,i^l 
d+l 

r+-= n 

i=0,*^l 



as d = An + 2 (46) 
as d = An (47) 



We note that in An + 2 dimensions for class BDI and 
CII and in An dimensions for class Dili and CI without 
enlarging the Dirac Hamiltonian, an extra mass term 



7i = < 



5]^7i asd = 4n + 2, 

i=0 
d 

iS TT 7j as d = An, 



(48) 



B. The rest of the classification 

The classification of the rests of the commutation and 
anticommutation relations has to be discussed case by 
case. The result of this classification for the four symme- 
try classes is shown in Table VI. 



can be found. It is to easy to check that 71 preserves 
system's non-spatial symmetries by Eq. (36). Further- 
more, 71 anticommutes with 7,^+1 but commutes with 
jd+i so that the reflection symmetry is preserved. This 
SPEMT 71 implies this system is classified by a '0' orig- 
inal topological invariant. 



1. Even spatial dimensions 

We consider systems in even dimensions. In the fol- 
lowing we show that the topological invariant in such a 
system is zero. For such a (i-dimensional system, in the 
Dirac Hamiltonian if only 70, 71, ... , 7<j, and S are in- 
gredients to construct reflection symmetry operator R, 
it is not possible that R satisfies the commutation and 
anticommutation relations in Table VI. To have R from 
the gamma matrices, we need to introduce an extra mass 
term 7 or an extra kinetic term "fd+i preserving the non- 
spatial symmetries. At the same time, the Dirac Hamil- 
tonian might be enlarged to have a new gamma matrix. 



2. Odd spatial dimensions 

In odd spatial dimensions without 7 and 7,2+1, we can 
use only the original gamma matrices to construct R sat- 
isfying the aforementioned commutation and anticommu- 
tation relations with T and C 

d 

=i ] | 7i, as d = An - 1, (49) 

i=0,i#l 
d 

= ] [ 7i, as d = An + 1. (50) 

1=0,1^1 
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The corresponding symmetry classes are class BDI and 
CII as d = An— I and class Dili and CI as d = 4n+l. The 
original classification of the AZ symmetry class shows 
that those classes have '0' topological invariants so that 
the non-spatial SPEMT 7 can be present in the mini- 
mal Dirac Hamiltonian without enlarging the dimension 
of the matrix. However, the anticommutation relation 
{7j li] = implies that 7 breaks the reflection symme- 
try. Such a system can be in a non-trivial phase. To find 
the topological invariant, we have to investigate the pres- 
ence of the SPEMT in the enlarged Dirac Hamiltonian 
in Eq. (14). 

All of the possible ways in Eq. (14) to doubling the 
size of the Dirac Hamiltonian can be separated into 
two expressions in Eqs. (20) and (21). To construct a 
SPEMT, all of the ingredients are the gamma matri- 
ces, a mass term 7, and the Pauli matrices. We find 
that the first summation over only odd number of the 
terms in Eqs. (20) and (21) the SPEMT 2Jt can be found. 
That is, D)l = (i) U°n%i 7«, ® On in Eq. (20) and DJl = 

(i)77i Iln^i 7«i ® a u m Eq. (21), where the Hermiticity 
is adjusted by the presence or absence of i and u = x or y 
to satisfy the non-spatial symmetries. The commutation 
and anticommutation [931, j ni ®I] = and {9ft, r y nj <8>I} = 
implies that 971 commutes with the enlarged reflection 
symmetry operator R' = (i) Y\ i=0 ^ 7» ® I- Therefore, 
the reflection symmetry is preserved. On the one hand, 
Y] n .-Ai k ni "f ni ® a z over an odd number of terms gives '0' 
topological invariant. On the other hand, an even num- 
ber provides a non-trivial topological phase. In short, 
this shows that Z topological property when k\ vanishes. 
Such a system has the 2MZ(MZ) 45 topological invariant. 

The cases that we have not discussed are class Dili and 
CI in An — 1 dimensions and class BDI and CII in An + 1 
dimensions. In the original classification without reflec- 
tion symmetry those symmetry classes possess Z and 2Z 
topological invariants. Let us first consider the Z case. 
To build the reflection symmetry operator R satisfying 
the sign changing of Rt and Rc in Eqs. (49) and (50), 
only using 70, 71, ... , 7<j, and S is not possible. There- 
fore, because o~ y plays a role in switching signs of Rt and 
Rc, we construct that the reflection symmetry operator 
is the direct product of a y and the reflection symmetry 
operator in the form of Eqs. (49) and (50) 

d 

R^ =i jj 7i<8>cr y , as d =An — 1, (51) 

i=0,i^l 
d 

R-+ = jf[ 7,0 0-^, asd=4n + l. (52) 

i=0,*^l 

At the same time, the Dirac Hamiltonian must be en- 
larged in this unique way 

d 

%l = m 7o ® I + k ili ® I (53) 



to preserve all of system's symmetries. From the prop- 
erties of the Z original topological invariant, a SPEMT 
is absent from this Hamiltonian. Therefore, to distin- 
guish this non-zero topological invariant, "Hf needs to be 
enlarged in the two possible forms 

=fcl7l Ol + y^ fenjn, ® I ® CT Z + ^ ^nj-frij ® ^4x4, 
rii n j 

(54) 

=^171 ® o z + fc »i7n, ® I ® o- z + ^ k n .-f nj <g> I 4x4 . 

rii nj 

(55) 

Again the original Z protects the non-trivial topologi- 
cal phases when an even number of a z in the Hamil- 
tonians. For the case of an odd number, although an 
extra mass term, which preserves the non-spatial sym- 
metries, is present, we have to confirm that the re- 
flection symmetry is also preserved so that the system 
is in the trivial phase. The expression of the mass 
term is 9JT = (i) Y\°^ 7 ni ® &x <g> a u for Eq. (54) and 

9tt = W7o Ilri^o 7»» o- u for Eq. (55), where the 

presence or absence of i keeps the Hermiticity of 9Jt and 
choosing u = x or y makes 9Jt preserve the non-spatial 
symmetries. It is easy to check that 9Jt commutes with 
^ = (*) IliLo M li ® °~y ® I- Hence, the behavior of 
even and odd numbers of o~ z shows that even with reflec- 
tion symmetry such a system still possesses Z properties. 
However, the reflection symmetry requires the dimension 
of the Dirac Hamiltonian to be doubled in Eq. (53) . The 
classification of the topological invariant changes from Z 
to 2Z. 

For the 2Z case, we know that in the proper basis the 
minimal Dirac Hamiltonian of 2Z is the two copies of the 
minimal Dirac Hamiltonian of Z: 

d 

R, 2Z = m 7 J ® I + J2 hjf ® 1 (56) 

The non-spatial symmetry operators can be expressed by 
the symmetry operator of Z 

Following the similar discussion of the Z system of 

the 2Z original topological invariant with the reflection 
symmetry is still classified as 2Z. 

VI. EXAMPLES 

Let us now discuss several examples of topological 
phases protected by reflection symmetry: more specifi- 
cally, we will consider: 

class AIII+i?+ and BDI+.R++ in d = 2: All of the 

non-spatial symmetry operators commutes with R. 
Therefore, Table III shows that both of the sym- 
metry classes possess MZ topological invariants. 
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class DUI+i? in d = 2: In this symmetry class, R 

anticommutes with T and C. This situation hence 
falls into Table IV. As described below, gapped 
phases in this case are characterized an MZ topo- 
logical invariant. 

class D+i? + in d = 2: With reflection symmetry, 
gapped phases in this symmetry class are charac- 
terized as a Z 2 topological invariant in Table III. 

class CII+i? in d = 2: For class CII no non-zero 

topological invariants are in the original classifi- 
cation or for mirror symmetry class A in the re- 
flection symmetric plane. The gapped phases are 
always trivial. 

class DUI+i? |_ in d = 3: R anticommutes with T but 

commutes with C. Hence, Table V shows there is 
a Z 1 in this case. 



A. Class AIII+i?+, BDI+.R++, and DUI+i? in 

d= 2 

Symmetry class AIII has two physical interpretations, 
one in terms of charged (complex) fermions with con- 
served fermionic number, and the other in terms of BdG 
Hamiltonians. As an electron system, a way to obtain 
such a system is to consider lattice fermion systems with 
bipartite hopping only. In this context, chiral symme- 
try of class AIII is sublattice symmetry. Alternatively, 
symmetry class AIII can be realized as a time-reversal 
symmetric BdG Hamiltonian with conserved S z spin ro- 
tation. While it is perhaps fair to say that the BdG inter- 
pretation is more experimentally realizable, since achiev- 
ing an exact sublattice symmetry is challenging, in this 
section we focus on electronic realizations of symmetry 
class AIII. 

Similarly to class AIII, symmetry class BDI has two 
physical interpretations, one in terms of charged (com- 
plex) fermions with conserved fermionic number, and the 
other in real (Majorana) fermions. Below, we will first 
discuss realization in terms of complex fermions; we will 
later discuss a realization in terms of Majorana fermions. 

Recall that in d = 2, there is no topological insulator 
in AIII and BDI if we do not impose reflection symme- 
try. Class Dili is of Z 2 type. With reflection, there are 
topological insulators in these three classes characterized 
by MZ topological invariants. 

a. bulk Hamiltonian Let us start by considering the 
following tight-binding Hamiltonian: 

ff = E^( r )(il !^)v>(r + a) + h.c. 

+ ^ (r) ( _* A A t ) V(r + V) + h.c. 

+ ^( r )(o -m)^ 7 ")' (58) 



where the two-component fermion annihilation opera- 
tor at site r, ip(r), is given in terms of the electron 
annihilation operators with spin up and down, c rl y 2 , 

as ip T (r) = (c Ti i, c rt 2), and we take t = A = 1 and 
/i = m + 2. The chiral p-wave superconductor has been 
discussed in the context of superconductivity in stron- 
tium ruthenate 46 and paired states in the fractional quan- 
tum Hall effect 47 . There are four phases separated by 
three quantum critical points at fi = 0, ±4, which are la- 
beled by the Chern number as Ch = > 4), Ch = — 1 
(-4 < /j, < 0), and Ch = +1 (0 < fi < +4). The non-zero 
Chern number implies the IQHE in the spin transport. 48 
In momentum space, 

H= ]T ^(k)[n(k)-a]i>(k), 

keBZ 

(-2Asinfc a; \ 
-2Asinfc !y . (59) 

2£(cos k x + cos k y ) + \i J 

A lattice model topological insulator in symmetry class 
AIII +R+, BDI +R++, and Dili +R can be con- 
structed by taking the two copies of the above two-band 
Chern insulator with opposite chiralities. Consider the 
Hamiltonian in momentum space, 

H = E E 4(k)[fi s (k)-a}Mk), (60) 
feeBZ s=t4 

where s =f, I represent "pseudo spin" degrees of freedom, 
and n s (k) is given, in terms of n as 

n-j-(fc) = n(k), n^(k) = n^(k) = n{k), (61) 

where k = (—ki, k 2 , ■ ■ ■) — {—ki,k±). I.e., 

H(k) = n x (k)r z cr x + n v {k)T a y + n z (k)T a z . (62) 

The model is chiral symmetric: 

S^nikjS = -H(k), S = t x g x . (63) 

The Hamiltonian is invariant under the following two 
TRS: 

T- 1 H(-k)T = H(k), T — t x o~qQ 7 T 2 = +l 
T~ 1 H(-k)T = H{k), T = T y a e, T 2 = — 1. (64) 

Also, the corresponding particle-hole symmetries with 
C = ST are 

C- 1 -H(-fc)C = --H(fc), C=T <j x e, C 2 = + l 

C- x U(-k)C = -7i(k), C=iT z a x <d, C 2 = + l. 

(65) 

Imposing the former form of TRS, the system falls into 
symmetry class BDI, whereas with the latter form of 
TRS, the system falls into symmetry class DHL 
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We now impose the following reflection symmetry: 

R- 1 U{k)R = 'H{k), R = t x . (66) 

Chiral and reflection symmetries commute with each 
other in the sense that 



[S, R] = 0. 



(67) 



For class BDI, R commutes with T and C (R ++ ). For 

class Dili, R anticommutes with T and C (R ). From 

Tables III and IV, both of the cases possess MZ topolog- 
ical invariants. Similarly, class AIII+i? + is also classified 
by a MZ topological invariant. 

b. MZ bulk topological invariant At the reflection 
symmetric plane, i.e., only k x — in the continuum 
model, H(k) commutes with R, and hence, it can be 
block diagonalized. Furthermore, since reflection R com- 
mutes with chiral symmetry in Eq. (67), each block has 
an off-diagonal structure in a proper basis, 



Tyro k ) ( n+ ^ k v) 
H ^' K y)-y o n-(o,k y ) 



D±{ky) 
Di(ky) 



(68) 



where ± indicates the eigenspace of R = ±1. For each 
block, the ID winding number (the topological invariant 
of symmetry class AIII in d = 1) is well-defined 18 ; by the 
definition of the MZ number in Eq. (27), we only need 
to focus on the winding number in one of the eigenspaces 
of i?, say +1, due to the absence of the weak index. The 
winding number is defined from q(k) in the Q-matrix 
in Eq. (A2). To construct the Q-matrix from the ( fey- 
dependent) occupied wavefunctions, we have to solve the 
eigenvalue problem 



(69) 



where a runs over occupied bands, a = 1,2 in our case 
and 



$3 = — , 



(70) 



Here we assume u a and v a are normalized so <E>± is also 
normalized. Because (H + ) 2 $± = (e a ) 2 $±, u a and v a 
are the eigenfunctions of D R D R and D R D R respectively 
and share the same positive eigenvalue (e a ) 2 : 



DD^u a = (s a ) 2 u a , D^D v a = (e a ) 2 v a . 



(71) 



Therefore, we can compute the projector of the occupied 
bands which have negative energies 



a ^ ' 

1(1 0\ _l v / u a v\ 
2 \0 l) 2^\v a ut 



(72) 



Therefore, due to Q(k y ) = I — 2P(k y ) and the definition 
of q(k y ) in Eq. (A2), we have the expression of 

q(ky)=J2\u a (k y ))(v a (k v )l (73) 

a 

the topological invariant for at k x = is defined by 

Vo = hj dky tr [« t ( fc «)^»«( fc «')] • ( 74 ) 

The winding number u n at the other symmetric point 
k x = 7r can be computed in the similar way so we have 
the MZ number N MZ by Eq. (27). 

Following the general discussion, we now calculate the 
topological invariant of the model in Eq. (60) . We first 
look for a basis where (i) R is diagonal, (ii) H(k) is block- 
diagonal for each reflection eigenvalue R = ±1, and (iii) 
chiral symmetry operation looks identical for both sectors 
of R. This can be achieved by the unitary transformation 
U = U1U2 where 



Ux = ^Oo^O +T z (T y ), U 2 = -^=(T +iT y ). 

Under the unitary transformation, 

URU^=r z a , USU^=r a z , 
UH(K)tf = ( %+{ n K) 



(75) 





n-(K) l ' 



where K = 0, 7r, and 



H%(ky) 





iriy =F n z 



-m y =F n z 




(76) 



(77) 



The topological invariant for the blocks satisfy 



4 



/'s —t/ K due ^° the absence of the weak index. In 



this specific model, i/q = 1 and 1/+ = as — 4 < \i < 
and i/q = and = 1 as < /U < 4. Therefore, by 
Eq. (27) N MZ = 1 and -1 as -4 < /U < and < fi < 4 
respectively. 

c. erfge theory Let us now introduce a boundary to 
the system. In the continuum model, one way to do this 
is to make the mass y-dependent m — > m(y). The edge 
Hamiltonian, in a suitable basis, is 



U{k x ) = k x a 3l {H{k x ), eri} = 0. 



(78) 



Or, in the second quantized language, the edge mode is 
described by 



H 



dx ip' L id x ip L - ip R id x ipf 



(79) 



Since there are left and right movers, one could give a 
mass to gap them out. We can write down such two 
masses, 



(80) 
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if we have not imposed any discrete symmetry. Let us 
now impose the two discrete symmetries. With chiral 
symmetry, the first mass will be eliminated: the first 
mass term can be written as m^a x ^, where <J/ := 
(ipL,ipFt) T - Since {m<T x ,cr x } ^ 0, this mass term is not 
compatible with chiral symmetry. On the other hand, 
the second mass is m^^cF y ^. Since {m 5 a y7 a x } = 0, 
this mass term is chiral symmetric, and allowed to exist. 
With reflection symmetry, we should be able to prohibit 
the second mass term. By definition, reflection should ex- 
change ipL and tpR, and observing reflection should com- 
mute with chiral symmetry, R is given by 



i> L {-x). 

(81) 



Observe that reflection and chiral symmetry commute, 



[R,a x }=0. 



(82) 



The first mass as well as the kinetic term is invariant un- 
der TZ. However, the second mass is not invariant under 
TZ. Therefore, with both chiral and reflection symmetries, 
the edge state is stable. 

We have treated the case with unit topological invari- 
ant. By increasing the number of edge channels with the 
same or different signs in Eq. (79), any mass term is still 
prohibited. We can consider cases with MZ topological 
invariants. 



B. Class ~D+R+ in d = 2 

For symmetry class D (i.e., generic BdG systems with- 
out any symmetry) in d = 1, the system is classified as a 
Z 2 topological superconductor. Therefore, a Z 2 topologi- 
cal superconductor can be realized in a symmetry class D 
system with reflection symmetry R ++ due to the upward 
shifting in Table III. As follows, we will describe this case 
in more details with an example. 

Observe that the pairing terms in x-direction for the 
spin up and spin down sector differ by sign. The model 
is invariant under reflection defined by 



TZc< 



r TZ~ x = c_ 



(s =t4). 



(83) 



Furthermore, for spin-one- half system, R 2 = — 1. With- 
out affecting the equation above, we performs a phase 
shift so that R = a x is hermitian. In particular, 
note that the pairing terms in rr-direction transform un- 
der reflection TZ as TZ: £ r A(c\ r c\ r+S; - c\ r+& c\ r 

-E,A(c 



c| r+£ c| r ), and hence TZ exchanges 



H l+w and H - P +i P - 

d. bulk spectrum With the periodic boundary condi- 
tion, we make use of the Fourier transforms which trans- 
form the Hamiltonian into 



0<k x <7r k y 

*L(M : = (4, fc . c i,-k, 4,*' c i-k)- ( 84 ) 



where the kernel % kx (k y ) is block diagonal in spin indices 
and given by 



^fex (ky) 
H-k x (k y ) 



n x (k y ) 

o ni(k y ) 



(85) 



where = 2t(cosk x + cosfc y ) — \i and = 

2A(— ssmk x +isink y ). The particle-hole symmetry op- 
erator C = t x Q commutes with R = a x . It corresponds 
to class D+i? ++ . At the reflection symmetric points 
k x = 0, 7r, following the general discussion, we take a 
basis in which R is diagonal: this is achieved by a uni- 
tary transformation, 



R -> URU~ 



U 



1 

71 



i i 

i -i 



Accordingly, the BdG Hamiltonian at k x 

/ Hk x =o,Tr(k y ), is transformed as 



(86) 

0,7T, 



Hk x (k v )^ 



^■k' ± (ky) 



where note H 



Ht(ky) 

o n^(ky) 
nUk y ) + nUky) 



(87) 



Ti], =0n (k y ) in this example. 
+1 and R = — 1 sectors are 



<-k x =0,w(ky) 
The Hamiltonian for R 

identical, and given by 

nj± (u \ — ( ^ k v Afe B \ 

£ kv =2t(±l + cosk y )-n, A ky =2iAsmk y , (88) 

where + for k x = and — for k x = it. 

e. bulk topological invariant Once we decompose the 
Hamiltonian into reflection symmetric/antisymmetric 
(R = ±) sectors, we can define a d = 1 topological in- 
variant for each R = ± sector, following Ref. 42: Let 
us consider class D systems in d = 1 with translation 
invariance. In momentum space, it can be written as 

a,/3 k y 

where a, j3 are a band index, and we are using a Majorana 
fermions \ a (k v ) to represent a class D system. Then the 
Z 2 topological invariant is given by 

(-1)" := sgn[Pf B(0)]sgn[Pf B(tt)] = ±1. (90) 

In our problem, this Kitaev Z 2 invariant can be defined 
for k x = 0, 7r and for each R = ± sector. For a given 
sector (R = +, say), by performing the transformation 
to real Majorana fermion operators 



A fc :=c l v +%,> 

K ■= (c- kk -cl k )/i, 



A fc := 



Afe, 
A', 



(91) 
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the Hamiltonian in Eq. (88) is given by 

tlk x =oA k v> - 4 2^ A -k y y £ ky 2AsinfcJ Afe «' 

(92) 

Let us focus on the Z 2 topological invariant at k x = 0. 
Thus, the fj-matrix at k y = 0, n is given by 

(93) 

and hence the topological invariant is 

(-iy° =sgn[PfBo(0)]sgn[PfB (7r)] 

= sgn(4i - fj,) ■ sgn(-/x) (94) 

Similarly, the topological invariant at fc^ = ir is 

(-1)"* = sgn(-4i - ft) ■ sgn(- M ) (95) 

Therefore, the Z 2 topological invariant in Eq. (28) for the 
entire reflection system is 

(-1)^2 = sgn (+4f _ M ) . sgn (_4f - (j,). (96) 

The presence of the topological nontrivial phase is in the 
region — At < \i < At. 

f. edge state A consequence of non-trivial topology 
in the bulk is the presence of edge modes. Let us now 
consider an edge of d — 2 class D topological supercon- 
ductor: 

S = J dxdr (VtHWf + . (97) 

where ip^ and ipi are a left-moving (right-moving) Ma- 
jorana fermion. A mass term imipfipi, is prohibited be- 
cause of a reflection symmetry: 

R : V>f -> "0-1, ^4. -> V'T- (98) 

In other words, with reflection the edge state is stable, 
and there is a bulk topological superconductor protected 
by reflection. 

Such protection by reflection symmetry is of Z 2 type 
as we can readily see by considering two copies of the 
edge theory: 

S = dxdT^2(ipti(-id x )i; t i + ipn(id x )ipn) . (99) 
•* i=i 

The corresponding Hamiltonian is 

H = j dxf T Ht>, H = -id x T (J z . (100) 

where ^> T = (^fi, ^4.1 > ^t2> ^4.2 ) T - Mass terms 
iip-fiipli and iip^ip^ are again prohibited. However, 
-2^1-1^4.2 + ^41^2) = * T Tj,cr a; *, and -2i(V'tiV't2 + 
^4.1^4.2) = i ^ T Ty(jQ'4' 1 are invariant under reflection and 
hence allowed to be added as a perturbation. 



C. Class CII+i?__ in d = 2 

In class CII a system is invariant under TS, PH with 
T 2 = C 2 = — 1. When we impose reflection symmetry 
into the system and require R to anticommute with T and 

C, Table IV shows that class CII+i? in 2 dimension is 

always in the trivial phase. In this section, we will pro- 
vide one example to show that the gapless edge states and 
the entanglement mid-gap states inevitably are gapped 
by the translational breaking term, which preserves all 
of the discrete symmetries, in Eq. (31). Consider the 
Hamiltonian in class CII+i? 

Ti = M70 + sinfci7i + sin/j 2 72, (101) 

where M — m + cos k\ + cos fc 2 , 70 = t z , 71 = T x a Xl and 
72 = T x Uyii x . The corresponding non-spatial symme- 
try operators are T = a y Q, C — T x fj, y Q. Physically, Tj, 
(Tj, and m can be treated as particle-hole, one-half spin, 
and pseudo-spin degrees of freedom respectively. With- 
out enlarging the matrix dimension, we have extra three 
kinetic terms 73 = T x a z , 74 = r y , and 75 = T x a y /i z . 
In addition, sin ^73,4,5 preserve TRS and PHS. There- 
fore, the reflection symmetry operator can be defined as 
R = *7i73 = fj/i which anticommutes with T and C. For 
such a reflection system, as — 2 < m < 2 the topological 
phase is non-trivial if the translational symmetry is pre- 
served. That is, the real spectrum gapless edge states still 
can be observed at the y direction edge and the entangle- 
ment mid-gap states are present when we make a spatial 
cut exactly in half at x = 0. However, those states are un- 
stable when the translational symmetry is broken. Con- 
sider the Hamiltonian % with 5 A = ^717374 = r y a y and 
a translational symmetry breaking term c91 in Eq. (31), 
where 91 = i7d+i7<i+27d+3 = T y a x ii z . As shown in Fig- 
ure 2(a) the edge states are gapped when c > 77. Further- 
more, Figure 2(b) shows that the mid-gap in the entan- 
glement spectrum is destroyed; this confirms that in the 
region — 2 < m < 2 the system is in the trivial phase. 



D. Class DUI+i? + in d = 3 

By reading Table V, in three spatial dimensions class 
DUI+i? |_ possesses Z 1 . In what follows, we demon- 
strate Z 1 topological properties by considering this re- 
flection symmetric system in class Dili in three dimen- 
sions. First, without reflection symmetry we introduce a 
three-dimensional time reversal invariant topological in- 
sulator possessing artificial chiral symmetry correspond- 
ing to class DHL We choose the Hamiltonian of such 
a system to be identical to the Bogoliubov-de Genncs 
Hamiltonian of 3 He-B phase 49 ' 50 . The system of 3 He-B 
is recently shown to be a topological supcrfluid preserv- 
ing TRS and PHS in three dimensions 15-17 . Now instead 
of PHS, we require our system to be invariant under chi- 
ral symmetry even if TRS is broken. Therefore, we write 
down the Hamiltonian of the 3D topological insulator 
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(a)Real spectrum: size L x X L y = 140 X 24 

Entanglement spectrum 
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(b)Size L x x L y = 140 x 8 



FIG. 2. (a) The energy spectrum for the y edge states near 
E = as a function of a translational symmetry breaking term 
(c), and (b) the entanglement spectrum around the entangle- 
ment eigenvalue 1/2 as a function of c. We consider the open 
boundary condition is the y discretion. The parameters are 
set to m = — 1, 5 = arcsin(0.l7r), and r\ = 0.057T = 0.157. As 
c > 0.157, the edge states start being gapped and the entan- 
glement mid-gap states begin to move out from 0.5. However, 
as c < 0.157, -Emm is nonzero and the entanglement eigenvalue 
is not 0.5 due to the finite size effect in the numerical stimu- 
lation. 



possessing TRS and chiral symmetry as 



Hr 








tp+ 


-tPz\ 
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-tPz 
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-e P 





\-tPz 


-tp+ 








vf/. 



(102) 



± ip y and a ap 



where = (a tp , a± p , b ip ), p ± = p x 
and b ap operators indicate the two different sublattices 
instead of particle-hole degree freedom in 3 He-B phase. 
Let t — 2 and in the lattice model the Hamiltonian can 
be written as 



"Hdiii = 



sin p x T x a z + sm p y Ty + sin p z T x a x 



(103) 



where e p — 



to — cos — cos p y — cosp z and Ti/tjj de- 
scribes the sublattice/spin degrees freedom. The expres- 
sions of the TRS operator and chiral symmetry operator 



are T = a y Q and S — r x a y Q so the pseudo-PHS op- 
erator is C — t. x Q. Since T 2 = — 1 and C 2 = 1, such 
a system belongs to class DHL Due to chiral symmetry 
preserving, by using the similar method in sec. VIAOb 
the winding number N% can be computed: 



-2 as |m| < 1 
V = I 1 as 1< | m | < 3 
elsewhere 



(104) 



In the following, we put our focus on to = 2. In the 
original classification, this winding number corresponds 
to the number of the gapless surface modes and the en- 
tanglement mid-gap states, which are intact against any 
symmetry preserving disorder. 

Now we introduce reflection symmetry, which changes 
x — > —x. There are two possible expressions of the 
reflection symmetry operators R. First, a ap and b ap 
exchange under reflection so R = r x a x (we require R 
to be hermitian). This reflection symmetry is broken 
when the Hamiltonian is in the expression of Eq. (103). 
Secondly, the sublattice is invariant under reflection so 
R = a x . The reflection symmetry is preserved in our 
system. Therefore, R anticommutes with T and com- 
mutes with C. (For 3 He-B, such relations change because 
R = t z <t x .) Since R commutes with C, the Hamiltonian 
at the reflection symmetry planes k x — 0, 7r corresponds 
to class D in 2 dimensions which has a Z topological in- 
variant. Therefore, the mirror Chcrn numbers can be 
computed in those symmetry planes. We perform a uni- 
tary transformation on the symmetry operators and the 
Hamiltonian 



<?x -> ov, a z -> -a x , e -> e 



(105) 



so that R — a z and the Hamiltonian is transformed to 



Dili 



H+ 
sin p x 



sin p x r x 
H- 



(106) 



where 'H± = e p r z + sinp y T y ± svnp z T x in the eigenspace 
of R = ±1. As to = 2, by Eq. (A6) the Chern numbers 
for H+ at k x = 0, ir are Vq = = respectively. 

Therefore, the mirror Chern number defined by Eq. (27) 
is 



N, 



A II. 



1 



(107) 



The topological phase of the system protected by the Z 
and MZ topological invariants is described by the Z 1 
number in Eq. (39). In this case, N%i = 1 correspond- 
ing to the number of the gapless surface states and the 
entanglement mid-gap. 

Now consider the two copies of the topological insula- 
tor Hamiltonian in a system but with a sign changing. 
The enlarged BdG Hamiltonian can be written as 



o/2 _ 
rt DIII — 



"Hdiii 









Hi 



Dili 



(108) 
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where we choose R'^m = R.diii(Px —> ~Px)- The two 
topological numbers for H^m are Nmz = 1 and N% = — 1 
and then for the whole system are Nm z = 2 and N% = 
so N z i = 2. Although the Z topological number vanishes, 
the topological phase is protected by Nmz- 

x Let T-L' m n in H^m become "Hdiii with p x — > 
—Px,Py —Py Although A^ MZ = for the entire sys- 
tem, Nz — 2 so A^i = 2. That is the reason that there 
are two robust gapless surface modes. 

Also we consider H' DIII = "Hdiii with only p y — > — p y . 
Both of the two topological numbers for H' DIII switch the 
signs. Therefore, the total numbers vanish: Nmz and 
N%i . This phase is trivial because we can find a SPEMT 

which anticommutes with R.~k m . 



VII. CONCLUSION 

Combining Tables III to VI, we write down the com- 
plete classification results (27 symmetry classes +R) 
in Table I. The classification of reflection symmetric 
topological insulators and superconductors still has the 
same spatial dimensional periodicities with the original 
Altland-Zirnbauer classification in Table II. The complex 
and real symmetry classes with reflection symmetry have 
the periodsof two and eight respectively. Although the re- 
flection classification tables seem to be complicated, the 
two ingredients can slightly simplify the classifications. 
To define these ingredients, we consider a d-dimensional 
system in an AZ symmetry class with reflection symme- 
try. We denote N d as a topological invariant of the orig- 
inal strong index without reflection symmetry as shown 
in Table II. Furthermore, we define N^ 1 as the (d— 1)- 
dimensional topological invariant in the corresponding 
mirror symmetry class (appendix B). By observing the 
reflection tables, some of the topological invariants of the 
strong index are determined by N d and N^ 1 as shown 
in Table VII. In other cases, there is an ambiguity to 
determine the strong topological invariants of reflection 
symmetric systems. Therefore, the approach of using 
minimal Dirac model provides a systematic way to find 
topological invariants. 

While topological insulators and superconductors pro- 
tected by a set of spatial discrete symmetries arc more 
fragile in general than those protected by non-spatial 
symmetries, they are still fairly relevant to realistic sys- 
tems. We list several realistic reflection symmetric sys- 
tems in Table I. 

For example, we, again, note that for the Z 2 TRS 
topological insulator in class All (without reflection), 

its Dirac representative has reflection symmetry R 

in three spatial dimensions. This is, from the point of 
view of TRS topological insulators, somewhat acciden- 
tal. However, as we discussed, this is related to the fact 
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TABLE VII. The topological invariants of reflection symmet- 
ric topological insulators and superconductors, and the strong 
topological index for d-dimensional topological states in the 
original periodic table I (Nd), and the d— 1-dimensional mir- 
ror topological invariant (iV d_1 ). In several cases, N d and 
TV* -1 still cannot completely determine the reflection topo- 
logical invariant. We have to go back to the minimal Dirac 
Hamiltonian method to determine the topological characters. 



that there is a topological distinction of ground states 
even without time-reversal when reflection symmetry is 
a good symmetry. The associated topological invariant 
is integer- valued(MZ), as opposed to Z 2 . 

In fact, many experimentally realized topological insu- 
lators such as Bi 2 Se3 "accidentally" have reflection sym- 
metry. From the discussion above, even breaking TRS, 
a surface Dirac cone will not go away if both the time- 
reversal symmetry breaking perturbation and the surface 
geometry respect reflection symmetry. For example, we 
consider a surface which respects the reflection symme- 
try. When an in-plane magnetic field, namely parallel 
to the surface, is applied to the system, the time rever- 
sal symmetry is broken. However, the reflection symme- 
try is still preserved when the direction of the in-plane 
magnetic field is tuned to coincide with the reflected di- 
rection. If this occurs, the surface Dirac cone is stable 
since it is protected by the preserved reflection symmetry 
even though the time reversal symmetry has been broken. 
These stable surface Dirac cone may be detected by STM, 
say. (ARPES may not be ideal if we use a magnetic field 
to break TRS). When the direction of applied magnetic 
field is away from the reflection direction, the system re- 
spects neither reflection nor time reversal symmetry. As 
a consequence, a gap should open in the surface states. 
We predict that the surface magneto-resistance with an 
in-plane magnetic field changes significantly, when the 
field direction is rotated and when the chemical poten- 
tial is close to the Dirac point. This is because the surface 
gap varies with the rotation of the field direction. 

We close with discussion on the effects of disorder 
on topological insulators and superconductors protected 
by reflection symmetry. In the ten-fold classification of 
topological insulators and superconductors, it has been 
proved useful to consider the boundary (edge, surface, 
etc.) Anderson localization problem: For a topological 
bulk, one should find a boundary mode which is com- 
pletely immune to disorder. In turn, once one finds such 
"Anderson delocalization" at the boundary, it means 
there is a topologically non-trivial bulk. Not only this 
bulk-boundary correspondence can be used to find and 
classify bulk topological phases in the absence of disor- 
der, it immediately tells us such topological phases are 
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stable against disorder; two phases which are topolog- 
ically distinct cannot be adiabatically connected by ei- 
ther spatially homogeneous or inhomogeneous deforma- 
tions. For topological phases protected by a set of spatial 
symmetries, stability against disorder is, in general, not 
trivial, since spatial inhomogeneity does not respect the 
spatial symmetries. One can still consider, however, sit- 
uations where the spatial symmetries are preserved on 
average. 31,51,52 Below, we will consider the stability of 
reflection-protected topological phases we identified ear- 
lier against disorder which is reflection symmetric on av- 
erage. 

Let us consider as an example the reflection symmet- 
ric topological insulator in symmetry class A in three 
dimensions (class A + R). (For other examples in two 
dimensions, see Ref. 31 .) For symmetry class A in three 
dimensions, we have a topological insulator protected by 
reflection symmetry, which is characterized by an inte- 
ger topological invariant. Let us consider the surface 
Hamiltonian: W{r) = Wo{r) + V(r) where r denotes the 
two-dimensional coordinates on the surface, Ho(r) is a 
kinetic term (the surface Dirac kinetic term). We have 
added a random perturbation V(r). The disorder- free 
part is reflection symmetric under R, which is reflection 
inherited from the bulk, R~ 1 Ho(f)R = Ho(r), while dis- 
order V(r) is not so. The reflection symmetry can be, 
however, imposed on average: R V(r)R = V(r), where 
7TT represents the quenched disorder averaging. We could 
approach this problem by means of effective field theories 
of Anderson localization, the non-linear sigma models 
(NLcrMs); they describe slowly varying degrees of free- 
dom in a disordered metal, which are related to a diffu- 
sion motion of electrons (called "diffusions" and "Cooper- 
ons" ) . When derived for the disordered surface problem, 
the action of the NLctM is given by 

Snl<tM = \ J d 2 r tr [d^Qd^Q] 

+ f cPre^tr [Qd^Qd v Q] , (110) 

where a matrix field Q(r) is the NLctM field Q e 
U(2AV)/U(AV) x U(N r ) and N r is the number of repli- 
cas; A the coupling constant of the NLuM, which is the 
strength of interactions among diffusons and Cooper- 
ons, and is inversely proportional the conductivity. The 
last term the action is the topological term (Pruiskcn 
term) ,which counts the non-trivial winding associated 
to 7r 2 [U(2Ag/U(AV) x \J(N r )} = Z. (Here we are con- 
sidering the real space which is topologically equivalent 
to a sphere.) In the absence of any discrete symme- 
try, O can take, in principle, any value, in which case, 
one can make electrons to be Anderson localized. While 
the action for generic value of O breaks reflection sym- 
metry (on the surface), = (integer) x tt turns out to 
be consistent with reflection symmetry. Moreover, when 
6 = (odd integer) x ir there is no Anderson localization. 
The NLctM (110) can be derived for the Dirac represen- 
tative of the surface mode, which consists of N flavors 



of two-component Dirac fermions when the bulk topo- 
logical invariant is A G Z. The theta angle is given by 
O = Nir. Therefore there is an even odd effect; when 
the bulk topological invariant is odd (even), the surface 
mode is stable (unstable) against disorder. This would 
mean that, in the presence of spatially inhomogeneity 
which nevertheless preserves reflection symmetry on av- 
erage, the topological distinction is not Z, but Z 2 . 
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While we were finalizing the paper, a preprint 53 ap- 
peared in which symmetry class D with reflection symme- 
try is discussed. One of their considerations corresponds 
to class D + R ++ in our general classification scheme. 

Appendix A: The calculation of Z number 

We review the calculation of the topological number 
for the symmetry classes possessing a Z or 2Z topolog- 
ical invariant. For the classification without reflection 
symmetry, in odd (d = 2n + 1) spatial dimensions such 
symmetry classes always have chiral symmetry, which is 
a key point to define the topological number. In even 
(2n + 2) spatial dimensions, the topological number is 
exactly the same with the n + 1-th Chern number. 

To seek the expression of the topological number in 
2n + 1 spatial dimensions, first we introduce the spectral 
projector onto the filled Bloch states and the "Q-matrix" 
(flat band) 40 by 

fW = Ek"W)KWI, Q(k) = i-2P(k). (Ai) 

d 

We note that P(k) 2 = P(k) so Q(k) 2 = I. If the sys- 
tem possesses the chiral symmetry, which is described 
by Eq. (4), with negative energy filled N + = The 
Q-matrix can be brought into block off-diagonal form, 

Q(k) = ( gt ° &) 9( Q fc) ) , q(k) € U(N_). (A2) 

in some basis. The topological number in 2n + 1 spatial 
dimensions is characterized by the winding number 

«*H-iM = fer^(i)" +1 / tr[(9- 1 *) 9n+1 ] 
(271+ 1)! 27T J BZ d=2 n +i 

(A3) 
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Here, BZ means the integration over ti-dimensional k- 
space. That is, the region is the first Brillouin zone in 
the lattice model. 

For the topological number in In + 2 spatial dimen- 
sions, first we define the non-Abelian Berry connection 
of the occupied bands 

A & \k) - Af{k)dk, = (u- & {k)\du7{k)) , (A4) 

where /i = 1, • • • , d, a, b = 1, • • • , N_. And the Berry 
curvature is defined by 

T &i (k) = dA hl + (A 2 f l (A5) 

The topological number is captured by the n+l-th Chern 
number 

= Ch„ +1 [.F] = 7 -^— f tr(^)" +1 

(A6) 

Appendix B: Symmetry classes in mirror planes 

In the mirror planes k\ — 0, tt, the Hamiltonian com- 
mutes with the reflection operator R 

[n kl=0 ,„,n\ = o. (bi) 

Furthermore, because R is hcrmitian and R = 1, the 
Hamiltonian can be decomposed to two diagonal blocks 
H± corresponding to the eigenspaces R = ±1 respec- 
tively. Now consider non-spatial symmetries, which de- 
termine the symmetry class of the system. If all of the 
non-spatial symmetry operators commute with R, the 
two block Hamiltonians H± belong to the same symme- 
try class. Furthermore, we name this symmetry class for 
W.± in the mirror planes as mirror symmetry class. How- 
ever, the mirror symmetry class changes when at least 
one of the non-spatial symmetry operators anticommutcs 
with R. The Hamiltonians R± are invariant under the 
i?-commuting symmetries but not invariant under the R- 
anticommuting symmetry. Therefore, the mirror symme- 
try class is determined by the i?-commuting symmetries. 
Table I shows the mirror symmetry classes for each possi- 
ble algebraic relation between non-spatial symmetry op- 
erators and R. 



Appendix C: The proof of the Z 1 number definition 

We will prove that the bulk topology of the Z 1 system 
is determined by the maximum value of \N%\ and \Nmi\ 
in general cases. We simplify the problem by considering 
the Dirac Hamiltonian with the coefficient of the mass 
term: 

m = M — (k x i <5i) 2 — fc 2 . (CI) 



In the proof of sec. V A 1 we only consider Si — 0. In gen- 
eral, Si can be any number in BZ. As Si — 0, 7T, Nz and 
Amz can be computed for such Dirac Hamiltonians. As 
Si is not at the mirror symmetry points, the Dirac Hamil- 
tonians only make a contribution to Nz- To achieve the 
proof of the Z 1 number definition, we consider the Dirac 
Hamiltonians in any possible distribution and determine 
the interplay between the phase of the non-trivial Dirac 
Hamiltonians, Nz, and N mi- 
Id. general, the MZ topological number is computed 
in the two symmetry planes k x — 0, n. To simplify the 
problem, we consider the weak mirror index in Eq. (25) 
vanishes so vqv^ < 0; hence, the MZ number is given by 
Nmi = v§ + v-n- However, with translational symmetry 
breaking by folding the BZ, the two symmetry planes 
collapse to the one symmetry plane k = 0. The MZ 
number of this d — 1-dimensional plane is v Q + v v . Since 
JVmz is invariant under density waves that connect these 
two points, the d — 1-dimensional Z number in the new 
k x = 7r plane vanishes. Thus, we still can discuss any 
MZ number only in the k x — plane without loss of 
generality. 

The Z topological number (Nz) for the entire system, 
in general, cannot be defined in the mirror symmetry 
planes. The discussion of the Z number can be sepa- 
rated to two parts: Si = 0, 7r and Si ^ or it. First, 
the part N% of the Z number is contributed from the 
Dirac Hamiltonians with Si = 0, 7r. The problem can 
be simplified by considering only Si — case. The rea- 
son is that the contribution of the Si = it can be moved 
to Si = by folding the BZ. Secondly, the other part 
N% of the Z number are contributed from Si or n. 
Due to the reflection symmetry, N% must be even. Out 
of the mirror symmetry planes, the bulk topology is los- 
ing the reflection symmetry protection. Therefore, if the 
Dirac Hamiltonians (±5j) possess Z numbers that differ 
by signs, the system combined by these two Hamiltonian 
is in the trivial phase. Hence, we put our focus on these 
Dirac Hamiltonians having the same sign of the Z num- 
bers. 

We define (+) as a pair of the Dirac Hamiltonian (±<5j) 
with an even number of the gamma matrices having the 
minus sign. In other words, the Z number in this case is 
2. Similarly, (— ) indicates an odd number of the minus- 
sign gamma matrices and N!^ = — 2. Consider a system 
possessing the Dirac Hamiltonians with (+, +) (see sec. 
VA1) and (-). Hence, N z = -1 and N M z = 1. The 
entire Hamiltonian for the system is written as 

H= H \ n ++ h l 1 b ki + U- (4 1+tf c fcl+4 + c f fei _ s c kl _ s )] , 

fel 

(C2) 

where H++ = m-^ + £\ =1 k^i and H- = J2 n , k ri z lm - 
X)n dd knjiny Due to Eq. (13) the Hamiltonian is in- 
variant under the reflection symmetry with the reflection 
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operator in the second quantization 

R = Y,( b l Rb ->^+ c l Rc -^- 



Instead of adding an extra (+,+), we add (— , — ) 



(C3) 



The SPEMT is found to prevent two of the three mini- 
mal Dirac Hamiltonians passing through quantum phase 
transition: 

91 + + 

m± = 7^ kl ° kl+s ± k > Ckl - s) + h ' c -' (C4) 

where 9T = (i) Yl^. 7«j an d we choose the presence or 
absence of "i" to preserve TRS and PHS. Moreover, in 
order to preserve the reflection symmetry +/— in 91 cor- 
responds to rij 7^ 1/rij = 1 in the product of 91 The 
entire Hamiltonian with 9T± and a real coupling coeffi- 
cient g is given by 

(U++ , 9 9T \ 
H± = *± 9<rf U- *± (C5) 
fcx V U-) 

in the basis of 

(C6) 

Hence, the first two Dirac Hamiltonian become trivial 
and the last Dirac Hamiltonian survives. The topological 
number is one, which equals to Nji in Eq. (39). 

Now we add an extra (+, +) into the original system. 
That is, Nz — and N m % = 2. The entire Hamiltonian 
with the SPEMTs is written as 



in the basis of 



o 

girt 
V o 



o 

H++ 
g 2 rt 




g x m 

,g 2 or o 

H- 

H- 



*± (C7) 



(a^ b kl -4=(c fel+(5 ± c kl s) -^{c kl+ sT c^-sfj 

(C8) 

The SPEMT coupling from the last Dirac Hamiltonian 
7~Lod& to the others is forbidden by the symmetries. By 
performing a proper unitary transition, only two Dirac 
Hamiltonian couple so that these subsystems become 
trivial. The two other Dirac Hamiltonian keep the origi- 
nal bulk topology structure. Hence, the topological num- 
ber is two, which is the maximum value of \Nj\ and 
\Nmj.\- 



U — = -TO70 - fci7i 



i=2 



kiji 



(C9) 



into the original system. In other words, consider the 
system with N% = Nmz = 0. The entire Hamiltonian 
with the SPEMT is given by 



^ ± = E $t ± 



(U— 9l 2\ 
H++ g 2 m 
g 2 rt U- 
\gi$ o H-J 



<$>4 



(CIO) 



where £ = i9T7o7i preserves all of the symmetries. All 
of the Dirac Hamiltonian couple so the entire system is 
in the trivial phase. As expected, the topological num- 
ber vanishes. Thus, only one (+, +) can couple with a 
pair Dirac Hamiltonian (— ) and then one Dirac Hamil- 
tonian survives and provides the non-trivial bulk topol- 
ogy. In general, consider A +j+ (+, +) Dirac Hamilto- 
nians and A_/2 pairs of (— ) Dirac Hamiltonians in a 
system. The number of the nontrivial Dirac Hamil- 
tonians is |A +!+ — A_/2|. Similarly, the number of 
the protected modes for the other cases separately are 
|A +; _-A+/2|, |JV_ t _— JV_/2|, and | A_ j+ -A + /2|, where 
N+/2 is the number of the pairs (+) Dirac Hamiltoni- 
ans. However, the discussion in sec. VAl (+,±) the 
number of the protected modes for the entire system is 
|A +!+ — N + - \ + |-/V_ i+ — JV_ s _ | . In this case, we replace 



N±,± A±, T - A_/2 



(Cll) 



Moreover, a system with the same number of (— ) and 
(+) is trivial. Thus, the total number of the protected 
modes is given by 

A+-A_, lAr A7 N+-N_. 



2 ' ■ « -.-r -•- 2 

(C12) 

which is the topological number for the Z 1 system. By 
applying those identities 

A°±A MZ = A±, + - A±,_, (C13) 
A z = N + — Nz = A z ° + A z , (C14) 

we can simplify the expression of the Z 1 number: 



N, 



+ + , ,A°-7V MZ + A; 



+ 



=Max(|^|,|A M z|)- 



(C15) 
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